Space Density of Optically-Selected Type 2 Quasars 



Reinabelle Reyes^, Nadia L. Zakamska^'^, Michael A. Strauss^, Joshua Green^, Juhan H. 
Krohk^, Yue Shen^, Gordon T. Richards^, Scott F. Anderson^, Donald P. Schneider^ 



O 

(N 



6 



> 



o 

oo 

o 



X 



ABSTRACT 

Type 2 quasars are luminous active galactic nuclei (AGN) whose central regions are obscured 
by large amounts of gas and dust. In this paper, we present a catalog of type 2 quasars from 
the Sloan Digital Sky Survey (SDSS), selected based on their optical emission lines. The catalog 
contains 887 objects with redshifts z < 0.83; this is six times larger than the previous version 
and is by far the largest sample of type 2 quasars in the literature. We derive the [OIII]5007 
luminosity function for < ipiii] < lO^^L© (corresponding to intrinsic luminosities up 

to M[2500A] ~ -28 mag or bolometric luminosities up to 4 x 10^'' erg s This luminosity 
function provides robust lower limits to the actual space density of obscured quasars, due to our 
selection criteria, the details of the spectroscopic target selection, as well as other effects. We 
derive the equivalent luminosity function for the complete sample of type 1 (unobscured) quasars 
and determine the ratio of type 2 to type 1 quasar number densities. Our data constrain this 
ratio to be at least ^ 1.5 : 1 for < LrniTTi < lO^-^L^ at z < 0.3, and at least ^ 1.2 : 1 for 



^[OIII] 



IO^Lq at 0.3 < z < 0.83. Type 2 quasars are at least as abundant as type 1 quasars 



in the relatively nearby Universe {z < 0.8) for the highest luminosities. 
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Active galactic nuclei (AGN) can be classified 
as type 1 (unobscured) or type 2 (obscured) based 
on the presence or absence of broad hydrogen 
and helium emission lines in their optical spec- 
tra. Unification models of AGN attribute the 
distinction between these two types to a differ- 
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ence in the observer's viewing angle to a nucleus 
surrounded by non-isotropic obscuring material 
(|Antonucci|[l993t lUrrv fc Padovanilll995l) . These 
models are well-established for low-luminosity, 
nearby AGN. However, their applicability to high- 
luminosity AGN (i.e., quasars, classically defined 
to be sources with bolometric luminosities greater 
than 10'*^ erg s~^) has long been controversial. 
Moreover, while it has long been known that ob- 
scured AGN dominate the l ow-luminosity popu- 



lation in the loc al Universe ([O sterbrock fc Shaw 
1988; Salzcr 198 9:lHuchra fc Bur g 1992; Ha o et al. 
2005bl : ISimpsonI 20051 ). the situation is less clear 



for high-luminosity AGN. This is in part because 
these objects are much rarer and more difficult 
to sample, since the AGN luminosity function de- 
creases steeply with luminosity. 

In this paper, we present a catalog of 887 
type 2 quas ars from the Sloa n Digital Sky Sur- 
vey (SDSS; lYork et al.l l2000l ). the largest sam- 
ple of type 2 quasars in the literature to date. 
They are selected based on their optical emis- 
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sion lines, have redshifts z < 0.83 and [OIIIJ5007 
emission line luminosities extending to ipiii] ^ 
1O"'^''L0 (corresponding to intrinsic UV luminos- 
ity M2500 ^ —28 mag). Multi- wavelength ob- 
servations of the most luminous objects from a 
previ ous version of the sample ( Zakamska et al.l 
2OO3L hereafter Paper I) have confirmed that they 
have infrared luminosities up to and above 10^^ 
erg s~^, have spectral en ergy distributions ex- 



pected of type 2 quasars ( Zakamska et al. 2004t 



Vignali et al 
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Ptak et al.ll200alVignah et al 



20061 : IZakamska et a l. 20081) , and contain type 1 
quasars in thei r centers revealed by polarimetric 



quasars m tnei r centers reveaiea by polari i 
measurements ( Zakamska et al. 2005 . 120061 ) 



With this updated catalog drawn from roughly 
three times as much SDSS data, we can now 
sample the high-luminosity AGN population in 
sufficiently large numbers to draw quantitative 
conclusions. First, we derive the [OIII]5007 lumi- 
nosity function (LF) of type 2 quasars. Then, by 
directly comparing the space densities of type 2 
and type 1 sources, we place robust lower lim- 
its on the fraction of obscured quasars in the 
local Universe as a function of [OIII]5007 lumi- 
nosity. Studying the space densities of differ- 
ent types of AGN provides strong constraints 
on the simplest AGN unification scenario as 
well as its modifications. Moreover, quantifying 
the obscured quasar population is essential for 
many applications, such as relating the present 
mass density of local black holes to the accre- 
tion history of the en tire AGN popul ation (e.g.. 



Sohan 1982; ,Yu fc Tren iainc 2002; Marconi et al 



20041 ). understanding the origin of the c osmic 
X-ray backgrou nd (e.g., IComastri et al. 19951 
Gilli et al.ll2007l ). and studying t he effects of lu- 
minosity on AGN structure (e.g., Lawrenc3 1991 



Urrv fc Padovanilll995l: iHopkins et al.ll2006l ) 



Our results are complementary to those de- 
rived from obsc ured quasars selec t ed from hard 
X-ray surveys jU eda et al. 2003t Szokolv et al 



2004t iBarger et al.. 2005; iMarkwardt et al. "1200a 



Treister et al .l2006HBeckmann et al.lbcioi] Sazonov et al 



2007) 



2005 



and mid -infrared c olor selection (|Lacv et al 
Stern et al,, .2005|: i Martmez-Sansigre et al 



2006[ [Polletta et al.l l2008l). For example, up to 
20% of hard X-ray selected AGN do not show any 



emission lines in their optical spectra (jRigbv et al 



20061 ) ■ and therefore would not be included in our 
sample. On the other hand, our sample would 



include Compton-thick objects that are missed 
in X-ray surveys, as long as their optical spectra 
meet our selection criteria. 

The paper is structured as follows. In Sj^l we 
discuss our selection of type 2 quasars based on 
their optical emission lines and in Sj3l we deter- 
mine the [OIII]5007 luminosity function from this 
sample. In SJH we determine the equivalent lumi- 
nosity function from a complete sample of type 
1 quasars and in Sj5l we determine the ratio of 
type 2 to type 1 quasars using the derived lumi- 
nosity functions. We discuss caveats and impli- 
cations of these results in ^J6] and we conclude in 
f|7l We adopt a 'concordance' cosmology, h = 0.7, 
r^M = 0.3, and 51 a = 0.7. We identify emis- 
sion lines using air wavelengths, identify objects 
with J2000 coordinates, and use asinh magni- 
tudes (ILupton et al. 19991) corre cted for Galactic 
extinction ( Schlegel et al.l 19981 ). We often ex- 



press luminosities in units of solar luminosities 
L© = 3.826 X 10^3 erg s'^. 

2. Type 2 Quasar Sample 

Our sample of type 2 quasars is selected from 
the SDSS spectroscopic database as objects with 
high-ionization, narrow emission lines, following 
Paper I. We describe SDSS data processing in 
the targeting of objects for spectroscopy in SDSS 
in H2.2\ our spectroscopic selection criteria in §2.31 
and measurement of the [OIIIJ5007 emission line 
in 



2.1. SDSS Data 

The SDSS has imaged - 10,000 deg^ of the 
sky with good astro n ietric a nd photometr i c cal- 
ibration (IPier et al.l 120031: Priacker et al.l l2006t 



1996 



Smith et al.l2002tllvezic et al.i2004!:'Padmanabhan et al 



2008 ) in the SDSS ugriz filters f|Fukugita et al 



Stoughton et aL 20021) using a dr ift-scanning 



wide-field camera ( Gunii et al.l 19981) on a ded 
icated 2.5 m telescope ( Gunn et al.l 20061) . For 
each object the photometric pipeline returns var- 
ious measures of fiux in each band, such as Pet- 
rosian (1976), poi nt spread function (PS F) and 
model magnitudes ( Stoughton et al. 20021) . 

A subset of objects from the imaging survey are 
targeted for spectroscopy and assig ned to a series 



of pla tes containing 640 fibers each (jBlanton et al 
20031 ). each of which subtends 3" on the sky; the 
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spectra cover 3800 — 9200A with resolution of 
1800 < R < 2100 and sampling of ~ 2.4 pixels 
per resolution element. The relative and absolute 
spectrophotometric calibration are good to ~ 5%. 
Spectral flux errors per pixel are typically of the 
order 1 x 10~^^ erg s~^ cm~^ A~^. 

Two independent spectral reduction pipelines 
assign redshifts and classifications to these spec- 
tra and measures fluxes of severa l major emission 



lines. The spectrold pip eline (|Stoughton et al 



20021 : ISubbarao et al.ll2002l ). fits Gaussian profiles 
to emission lines to determine emission line red- 
shifts, and the specBS pipeline, written by D. 
Schlegel, carries out fits of spectra to templates 
in wavelength space. For 98% of spectra, the 
measured redshifts from the two pipelines agree 
within 300 km s~^ for galaxies, and 3000 km 



for quasars ([Adelman-McCarthv et al.ll2007^ : but 
some type 2 quasars are among the 2% of the dis- 
crepant objects and we must treat them separately 



2.2. SDSS Spectroscopic Target Selection 

Spectroscopic target selection in the SDSS is 
based on a combination of photometric proper- 
ties, such as magnitudes, colors, and morpholo- 
gies, and in some cases, radio and X-ray prop- 
erties. In this section, we describe the subset of 
target algorithms from the 'Main survey' and the 
'Special Southern Survey' which are important in 
the selection of type 2 quasars, and subsequently, 
for our calculation of their luminosity function in 

SI 

The Main survey constitutes close to com- 
plete samples of galaxies and quasars, which com- 
prise roughly 75% of all SDSS spectra. The 
Galaxy algorithm targets resolved sources down 
to a limiting Petro sian magnitude of r = 17.77 
(jStrauss et al.l 120021 ) after correction for Galactic 
extinction following ISchlegel et al. I (|l998l) . These 
objects represent about 62% of the objects in 
the Main survey. The Luminous Red Galaxy 
(LRG) algorithm targets sources based on the 
distinctive colors of LRG s, down to a Petrosia. n 



magnitude of r = 19.5 (jEisenstein et al.l 120011) . 



representing 9% of the Main survey. The Low-z 
QSO algorithm targets mostly UV excess sources 
down to ipsF — 19.1, and the High-z QSO al- 
gorithm targets point sources with the colors of 
high-redshift (z > 3) quasars down to ipgp = 20.2 



([Richards et al.l 120021 ) . Taken together, quasar 
candidates represent 13% of the Main survey. 

Other Main survey target algorithms are as- 
signed spectroscopic fibers only after the galaxy 
and quasar targets have been allocated, so they 
do not produce complete samples. Objects with 
unusual colors, radio emission det ected by the 
FIRS T survey ('Serendipity FIRST'; 'Becker et ajj 
19951 ). and X-ray em ission detected by t he ROSAT 
survev f'ROSAT': IVoges et al.l Il999l see also 
Anderson et al. 20071 ) are selected down to a fibe r 
magnitude oi i = 20.5 (|Stoughton et al.l 120021 ) . 
These so-called 'serendipity' targets constitute 
about 5% of the Main survey. 

Another 10% of spectra are taken in the Equa- 
torial Stripe of the Southern Galactic Cap, cover- 
ing ~ 3 00 deg^, as part of the Special Southern 



survey (jAdelman-McCarthy et al.l l2006l ). Spec 



troscopic targeting algorithms for this survey are 
somewhat different from those in the Main sur- 
vey. Those important for the selection of type 2 
quasars are: (i) modified versions of the Main sur- 
vey algorithms, which target galaxies and quasars 
to limiting magnitudes < 1 mag fainter, (ii) the 
'faint quasars' algorithm, which is a modified ver- 
sion of the Main Low-z QSO algorithm with looser 
color criteria, and (iii) the 'photoz' algorithm, 
which targets sources in the blue end of the nor- 
mal galaxy distribution down to a Petrosian r- 
band magnitude of 19.5. These algorithms are 
exploratory in nature and have changed several 
times over the course of the survey. For this rea- 
son, we do not use objects targeted by these "in- 
complete" algorithms in the calculation of the lu- 
minosity function. 

2.3. Spectroscopic Selection of Type 2 
Quasars 

Following Paper I, we select Type 2 quasars 
as objects with narrow emission lines without un- 
derlying broad components, and with line ratios 
characteristic of non-stellar ionizing radiation. We 
searched the entire SDSS spectroscopic database 
as of 2006 July, which contains 1.08 xlO^ spectra 
from 1770 plates, before accounting for duplicates. 
This correspond s to ^ 80% of the DR6 spectro - 
scopic database (jAdelman- McCarthy et al.ll2007l ). 



and is three times larger than the one used in Pa- 
per I. 
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The type 2 quasar selection was performed us- 
ing the spectrophotomet ric cahbration used i n 



the First Data Release (lAbazaiian et al.l l2003r) 



The Second Data Release paper (jAbazaiian et al 



20041 ) describes a substantial improvement in the 
spectrophotometric calibration algorithm, which 
affects the ratios of line strengths between the red 
and the blue end of the spectrum. An additional 
modification in the algo rithm was incorporated in 



the Sixth Data Release ( Adelman-McCarthv et al 



20071 ) ■ whereby spectra are calibrated relative to 
PSF rather than fiber magnitudes of the standard 
stars measured on each plate. Thus, the flux scale 
used in our luminosity calculations f §2.4l and 
is on average 38% higher than what was used in 
the initial selection. These changes in calibration 
cause some minor incompleteness at the faint end 
of the luminosity function, but very few objects 
change their classification from AGN to star form- 
ing galaxy (Eqs. 1-3 below) due to these changes. 

Our automated selection algorithm applies a 
series of constraints on the data: a redshift max- 
imum, a spectroscopic signal-to-noise ratio (S/N) 
minimum, an [OIII]5007 luminosity minimum, 
and a set of emission line ratio cuts. The resulting 
~ 4000 spectra are then visually inspected and 
fit for broad components in H/3 and Ha for final 
selection. Redshifts are also checked for accuracy. 

The major differences in our selection algorithm 
from that in Paper I are: (i) we now impose a 
luminosity cut, i[oiii] > lO^ '^i©, because selec- 
tion based on emission line ratios becomes incom- 
plete and noisy at low luminosities; (ii) we consider 
all objects with redshifts z < 0.83 (our previous 
sample was at redshifts 0.3 < z < 0.8 3, comple- 
mente d b y the z < 0.3 AGN samp les bv lHao et al. 
2005al and lKauffmann et aL 1 120031) : and (iii) we use 



an improved algorithm for identifying weak broad 
components in Ha and H/3 (Fig. [Hand text below). 

We restrict the selection to objects with red- 
shifts z < 0.83 so that the [OHI]5007 line (the 
strongest expected emission line) is present in 
all spectra. In order to select emission line ob- 
jects, we require the rest-frame equivalent width of 
[OHI]5007 to be > 4A. In addition, the signal-to- 
noise ratio must be > 7.5, where the signal is the 
flux density in the seventh brightest pixel over the 
entire spectroscopic range (3800-9200 A; about 
3840 pixels) and the noise is the median estimated 
flux error per pixel over all pixels. This unconven- 



tional criterion allows retention of objects with 
weak continua but strong narrow emission lines, 
while rejecting continuum-dominated sources with 
low signal-to-noise ratio (Paper I) . 

When the redshift and classification of the two 
reduction pipelines spectrold and specBS (de- 
scribed in i 32.1|) agree, we apply different emission 
line criteria depending on the redshift of the ob- 
ject. For objects with redshifts z < 0.36, both 
H/3-|-[OIII]4959,5007 and Ha-h[NII] 6548,6583 line 
complexes are covered by the spectroscopic data, 
and the classical emission line diagnostic diagrams 
( Baldwin et"al]ll98ltlVeilleux fc Osterbrocklll98"7 



Osterbrock 19891 ) are used to distinguish between 
a stellar and an AGN ionizing continuum. The 
Mgll 2799 emission doublet is either not covered 
by the spectrum or falls into the UV where the 
spectra typically have low signal-to-noise ratio, 
so this line is not used for this set of spectra. 
We use l ine diagnostic criteria of the form sug- 
gested bv lKewlev et al. Idlool,) to distinguish type 
2 quasars from star-forming galaxies and narrow- 
line AGN. We require the ratio of luminosities 
7^ = [OIII]5007/H^ to satisfy either: 



log (7^) > 



0.61 



log([NII]6583/Ha) -0.47 



1.19 (1) 



or 



log(7^) > 



0.72 



log([SII]/Ha) - 0.32 



+ 1.30, (2) 



where [SII] refers to the combined luminosity of 
the doublet [SII]6716,6730. 

For objects with redshifts 0.36 < z < 0.83, the 
Ha-|-[NII] line complex is not covered by the SDSS 
spectra, so the classical diagnostic diagrams can- 
not be used. We adopt the [OIII]5007/H/3 ratio 
requirement in reduced form: 

log (7^) > 0.3, if H^ is detected with S/N > 3 

(3) 

or we require that H/3 is undetected, while 
[OIII]5007 is detected. In addition, for z > 0.6, 
the full width at half maximum (FWHM) of Mgll 
2799 is required to be < 2000 km/s. 

In the 2% of cases for which the two spectro- 
scopic reduction pipelines do not agree on redshift 
or classification, we find the emission line closest 
to the expected position of [GUI] 5007, given ei- 
ther of the two redshifts. If the line satisfies the 
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appropriate equivalent width and luminosity cri- 
teria above, we retain the object as a candidate. 

These selection criteria are designed to be max- 
imally inclusive. In particular, the line diagnostic 
criteria given in Eqs. [1] and [2] are applied with 
the 'OR' operator, in case some of the lines are 
not measured properly. At redshifts z > 0.36, 
only a very mild line ratio cut (Eq. [3]) is im- 
posed. Apart from the Mgll 2799 width crite- 
rion on z > 0.6 objects, no criteria to explic- 
itly reject broad-line AGN are imposed at this 
selection stage, since the weak H/3 lines com- 
mon in type 2 AGN are often poorly measured 
by the spectroscopic pipelines. Nevertheless, the 
vast majority of type 1 AGN are easily rejected 
by the automated procedure; because of their 
broad Ha and H/3, they tend to have small val- 
ues of log([OIII]5007/H/3), log([NII]6583/Ha) and 
log([SII]6716, 6730/Hq!), making these objects fail 
the line diagnostic criteria given in Eqs. [T] and [21 
This automated procedure selected ~ 4,000 ob- 
jects for visual inspection. 

At the visual inspection stage, we pursue two 
goals: (i) to remove objects with broad com- 
ponents in their permitted emission lines (the 
major contaminant at z < 0.36); (ii) to reject 
star-forming galaxies (the major contaminant at 
z > 0.36, since at z < 0.36 the diagnostic dia- 
grams are robust). We consider the [NeV] 3346, 
3426 emission lines to be una mbiguous signs of an 
underlying AGN continuum ( Groves et al.l 12004 ; 
Nagao et al. I l2006l) . so if either one of these lines is 
detected (if a single line is detected, it is usually 
[NeV] 3426), the object is considered an AGN. If 
neither of the [NeV] lines is detected, we use the 
criterion FWHM aOIII]5007)>400 km/s as a sign 
of AG N activity ( Zakamska et al. 2003t Hao et al 
2005a[ ). Since the [NeV] lines are weak, our se- 



lection is more robust at high luminosities, where 
these lines are more likely to be detected. 

Objects identified as AGN are then checked for 
the presence of broad lines. A broad component 
in the H/3 line can be rejected by examining the 
difference in of a single Gaussian and a dou 



ble Ga ussian fit to this line, similar to iHao et al 



(|2005aD . However, identifying weak broad compo- 
nents in Ha is complicated by the non-Gaussianity 
of emission lines and blending with [Nil] 6548, 
6583. Thus, we do not automatically reject objects 
for which a four-Gaussian fit is statistically pre- 



ferred to a three- Gaussian fit to the Ha-I-[NII] line 
complex. Instead, we use the non-parametric line 
fitting procedure iUustrated in Fig.[T] First, we de- 
rive narrow line profiles from the [OIII]4959, 5007 
emission lines, and then fit the Ha-|-[NII] complex 
assuming that all three lines have the same pro- 
file (with the ratio [Nil] 6583/ [Nil] 6548 fixed to 3). 
We visually examine the residuals, and if the non- 
parametric fit can reproduce the width of the com- 
plex, we retain the object in our sample. 

This selection procedure results in the sample 
of 887 type 2 quasars given in Table [H Figure [2] 
shows sample spectra of high [OIII]5007 luminos- 
ity objects in our sample, in order of increasing 
redshift. More than 90% of the type 2 quasar can- 
didates that were selected in Paper I are recovered 
by the automated selection procedure described 
above, so we estimate that roughly the same per- 
centage of type 2 quasars in the SDSS database 
are successfully selected. The spectra of those ob- 
jects from Paper I that we did not recover tend to 
be of low signal-to-noise ratio or have ambiguous 
classification. 

2.4. Measurement of [QUI] 5007 Luminosi- 
ties 

We measure [OIII]5007 line luminosities from 
spectra calibrated with the most recent spec- 



troph otometric calibration algorithm (jAdelman-McCarthv et al, 
20071 . also see beginning of We do not 

correct the measured line luminosities for dust 
extinction because of the large uncertainties in- 
volved in such a correction; we discuss how this 
effect might affect our results in We fit the 

[OIII]4959,5007 lines with two Gaussians of the 
same width and a fixed 1 : 3 amplitude ratio, plus 
a linear continuum, over the wavelength range 
4860-5060A. We note that Fell emission tends to 
be negligible in type 2 quasars (unlike in type 1 
quasars, for which we carefully subtract Fell emis- 
sion before measuring [OIII]5007 luminosities, as 
described in i j4.2l) . 

To check the fits, we also obtain a non- 
parametric measure of the [OIII]5007 line lumi- 
nosity by integrating the detected flux density and 
subtracting the continuum contribution. As Fig. [T] 
shows, this line is often significantly non-Gaussian, 
especially at its base. Asymmetries in line profiles, 
thought to be caused by radial outflows in the 
narrow-line region, have been observed in many 
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AGN (|Heckman et all [ToSlt IWhittld Il985l) . We 
find that the non-parametric measure of the lumi- 
nosity is systematicahy larger, by 5% on average, 
than the Gaussian measure; this indicates the 
presence of non-Gaussian wings in the [OIII]5007 
profiles of objects in our sample. 

Both Gaussian and non-parametric measures 
are given in Table [TJ In the rest of the paper, 
iyjojji] refers to luminosities measured from Gaus- 
sian fits. About 30% of the sample (257 objects) 
have i[oiii] > lO^'^L©. As we explained in 
the luminosity cut we imposed in our selection, 
-^[oiii] > 10^ '^ -^07 'w^s applied to luminosities cal- 
ibrated using a fiux scale that is ^ 38% lower than 
the one we use for these line measurements. Due 
to this change in spectrophotometric calibration, 
there are 149 objects with measured [QUI] 5007 
luminosities below of these, around 80% 

are within 30% of this value. 

3. Type 2 Quasar Luminosity Function 

In this section, we derive the [OIII]5007 lumi- 
nosity function from our sample of type 2 quasars. 
Since type 2 quasars are heavily obscured, the 
optical continuum luminosity is by definition a 
poor indicator of their bolometric luminosity. The 
[OIII]5007 emission line offers a plausible alter- 
native, since it is clearly detected in both type 1 
and type 2 AGN, and is observed to correlate with 
continuum luminosity for type 1 sou r ces (see dis- 
cussion in MA\ and ^6.1^ . iHao et aL I (l2005bll and 
SimpsonI (|2005l ) derived the [OIII]5007 luminosity 



function from a sample of type 2 AGN with red- 
shifts z < 0.3 and [OIII]5007 luminosities up to 
~ 10^-^ Lq and 10^-^ Lq, respectively. Our sample 
probes both larger redshifts and higher luminosi- 
ties, up to z = 0.83 and i[oiii] ~ lO^^L©. To 
minimize the systematic effect of redshift evolu- 
tion within the sample, we present the LF for three 
ranges in redshift: z < 0.30, 0.30 < z < 0.50 and 
0.50 <z < 0.83. 

The type 2 quasar LF we derive is a lower 
limit to the true LF for several reasons. First, 
obscured quasars for which the narrow emission- 
line region is also obscured would not meet our 
selection criteria, so they a re not included in our 
sample ( Rigbv et al. 20061) . Second, there may 
be objects that fall outside the regions in the 
multi-dimensional parameter space of colors, mag- 



nitudes and morphologies, which the SDSS spec- 
troscopic target algorithms are designed to target, 
so they are also not included in our sample. Third, 
we use line luminosities that are not corrected for 
dust extinction as it is unclear how to do so con- 
sistently (see %.2\i . Correction for line extinction 
would shift the LF toward higher luminosities and 
would yield a higher space density at any given 
luminosity. 

In H'S.1\ and ii3.21 we discuss the calculation of 
the LF, and in ^3.3[ we present the results and 
compare it with previous work. The derived type 
2 quasar LF is shown in Fig.[5l In §3.41 we discuss 
the efficiency of our selection of type 2 quasars and 
its implications to the interpretation of the derived 
LF. 

3.1. 1/Vinax Luminosity Function 

To derive the lu minosity funct ion, we use the 
l/Knax estimator ( Schmidtl Il968[) . in which the 
contribution of each object to the luminosity func- 
tion is weighted by its available volume Knax- To 
calculate Knax, we need to determine how each ob- 
ject had been selected for spectroscopy. The SDSS 
spectroscopic target selection algorithms that are 
important in selecting type 2 quasars are described 
in g22I They are listed in Table [2] together 
with the number of objects targeted by each al- 
gorithm. An object can be targeted by multi- 
ple target algorithms, so these numbers do not 
sum to the total number of objects in the cata- 
log. As noted in §2.21 target algorithms of the 
Special Southern survey have c hanged over time 
(jAdelman- McCarthy et al.ll2006l ). so we do not in- 
clude objects from this survey in the calculation 
of the LF. Furthermore, we include only objects 
targeted by any of the four primary Main sur- 
vey algorithms (Galaxy, Low-z QSO, High-z QSO 
and Serendipity FIRST) , which together represent 
83% of the fuU catalog (740 out of 887 objects). 
Fig. [3] shows the distribution of [OIII]5007 lumi- 
nosities and redshifts of objects targeted by these 
primary algorithms. 

In calculating Vmax, the line luminosity and 
S/N criteria we have applied in W2.'3\ are not im- 
portant for the luminosity and redshift ranges of 
interest here, since the [QUI] 5007 line is strongly 
detected in all objects in our sample. We there- 
fore focus on the selection criteria from each spec- 
troscopic target selection algorithm. For our pur- 
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poses, T4iax is given by the comoving volume over 
which an object would be selected by any of the 
four primary target algorithms: 

47r L, dz 



(4) 



where is the effective survey area, Zi and Z2 are 
the edges of the redshift range for which we are 
calculating the luminosity function, {dVd dz)dz is 
the comoving volume element in the redshift inter- 
val dz, and Q{z) is a generalized step function that 
is equal to zero if the object is not selected by any 
of the four algorithms we consider, and nonzero 
otherwise. We describe in detail how we calculate 
e(z) in go 

The effective survey area $7 is given by the area 
covered by the Main survey spectroscopic plates 
from which we have selected objects. To calcu- 
late ri, we define areas of intersection of plates 
and tiling rectangles and take th e area covered by 
their union (see Appendix B of IShen et al. 2007 
and references therein for details). This calcula- 
tion yields ~ 6293 deg^. Since some post-DR5 
plates are included in our parent sample, this area 
is slightly bigger than t he DR5 spectroscopic foot 



print area of 5740 deg^ (jAdelman-McCarthv et al 
20071) . 



The l/Knax luminosity function and its uncer- 
tainty are given by 



- 1 / 1 \ 



(5) 



and 



am 



1 \ ^^iax,7 



1/2 



, (6) 



where Lk is the mean [OIIIJ5007 luminosity of ob- 
jects in the fcth luminosity bin {L, L + {AL)k), and 
the sum is over the Nk objects in that bin. Fol- 
lowing common practice, we present our results 
in terms of the number of quasars per unit vol- 
ume per logarithmic luminosity interval, denoted 
byl>(L) = (L/logioe)ci>(L). 

3.2. Calculation of the Selection Function 

e(z) 

In this section, we calculate the function Q{z), 
which appears in the expression for VJnax given 



above (Eq. |4]). We determine whether an object, 
with redshift Zobs, would be selected by any of the 
four primary target algorithms, if it were placed 
at some other redshift z, based on its intrinsic 
spectral energy distribution (SED). For this ob- 
ject, Q{z) is given by a generalized step function: 



e(z) 



if selected by any of G, L, or H 
if selected solely by S 
if not selected by G, L, H, or S 

(7) 

where G, L, H, and S stand for the Galaxy, Low-z 
QSO, High-z QSO, and Serendipity FIRST algo- 
rithms, respectively, and a is a numerical factor 
(smaller than 1) that weights the contribution of 
the radio-selected objects to the LF. We discuss 
the determination of a in i i3.2.4l In practice, we 
determine 0(2) in steps of redshift, from z = 0.1 
to 0.8, separated by Az = 0.01, and we approx- 
imate the integral in Eq. [4] as a sum over these 
redshift slices. 

To determine &{z) for a given object, we calcu- 
late what its observed optical and/or radio prop- 
erties would be if it were placed at redshift z, and 
apply the relevant selection criteria for each al- 
gorithm. For a ll objects, we apply k-corrections 
( Sandagel 196l[ l based on their observed SDSS 
spectra, which take into account both the shape 
of the SED and the shifting of emission lines in 
and out of bandpasses. The latter effect is impor- 
tant because objects in our sample have equiva- 
lent widths as large as 1,400 A (Paper I). To cal- 
culate k-corrections for the SDSS u and z bands, 
we extrapolate the observed spectrum outside ~ 
3,000—10,000 A using a constant flux density; our 
results are not sensitive to this extrapolation. In 
the following subsections, we describe details in 
the procedure for determining Qiz) that are spe- 
cific to each target algorithm. 

3.2.1. Main Galaxy Target Algorithm 

The Galaxy algorithm targets objects that have 
Petrosian magnitude r < 17.77 mag, r band Pet- 
rosian half-light surface brightness /150 < 24.5 mag 
arcsec~^, and satisfy a star /galaxy separation cut 



and a fiber magnitude cut ( Strauss et al ] |2002[ ). 

Since surface brightness changes only by a fac- 
tor of ~ 3 throughout the redshift range probed 
by the Galaxy algorithm [z < 0.3; Fig. [3]), and no 
object in the sample is close to the star/galaxy 
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separation cut, the magnitude cut alone deter- 
mines whether an object is selected. We consider 
an object to be selected by this algorithm if its 
scaled and k-corrected Petrosian r-band magni- 
tude is brighter than 17.77 mag. We assume for 
simplicity that the spectral shape within the Pet- 
rosian aperture is the same as within the spectro- 
scopic fiber. 

3.2.2. Low-z QSO and High-z QSO Target Algo- 
rithms 

Low-z QSO targets must have 15.0 < ipsF < 
19.1 and satisfy various color criteria, i.e., they 
must be outliers from the ugri stellar locus, have 
non-galaxy colors if they are extended, and not oc- 
cupy any of several exclusion regions in color space 
designed to eliminate whit e dwarfs, A-type stars . 



and other contaminants (jRichards et al.l 120021 ). 
All UV excess sources {u — g < 0.6) that are not 
in the white dwarf exclusion region and satisfy the 
magnitude limits are also targeted. 

The High-z QSO target algorithm is designed to 
recover quasars at redshifts beyond ~ 3.0, as the 
Lyman break moves across the photometric bands 
with increasing redshift. High-z QSO targets must 
be point-like sources with 15.0 < ipsF < 20.2 



(jRichards et al.ll2002[) . They must also be outliers 



in the griz stellar locus, or occupy certain regions 
in the color space where high-redshift quasars are 
expected to lie. 

To determine whether an object would be se- 
lected by these algorithms at a redshift z, we run 
the scaled and k-corrected PSF magnitudes and 
errors in all bands through the fin al version of the 



QSO target selection code v3_l_0 ([Richards et al 
20021) . 



3.2.3. Serendipity FIRST Target Algorithm 

Serendipity FIRST targets are sources that 
have fiber magnitudes 14.0 < g,r,i < 20.5 and 
have counterparts in the FIRST catalog within 
2" of the optical position. Radio sources are in- 
cluded in the FIRST survey catalog if they have 
peak 20 cm radio flux density S'20cin > 5a + 0.25 
mJy/beam and iS'20cm > 1 mJy/be am, where a 



is the local RMS noise in the field ([Becker et al 
Il995[ ). 

To calculate the fiber magnitude as a function 
of redshift, we use the azimuthally-averaged radial 



su rface brightness pro file (following Appendix B 
of IStrauss et al.ll2002[ ) to calculate the light that 
would fall into the aperture at each redshift, ne- 
glecting the effect of seeing. To calculate the ra- 
dio flux density as a function of redshift, we as- 
sume a power-law r adio spectrum Fy = Av ^ with 
^ = -0.5, following IZakamska et aD ([2004 ). Most 
of our objects are point radio sources, so we ap- 
proximate the redshift scaling of the peak flux 
density 5200111 to be the same as the total flux 
F^ cx Dl^{z){l + zf+f', where Di^{z) is the lu- 
minosity distance at redshift z. 

Since most objects targeted by the Serendipity 
FIRST algorithm were selected in the j-band and 
because the volume corresponding to the bright 
limit is negligibly small, we consider an object to 
be selected by this algorithm if it satisfies both the 
i-band fiber magnitude faint limit and the radio 
fiux limits. 



3.2.. 



Probability of Radio Detection 



About 30% of the objects included in our LF 
calculation are selected solely by the Serendip- 
ity FIRST algorithm. These objects do not have 
the optical morphologies or colors that would al- 
low them to be selected by the three other algo- 
rithms we consider. For this set of objects, we 
consider two definitions for the selection function 
@{z) — a(z,L[oni]) depending on whether an ex- 
plicit correction for the probability of radio detec- 
tion pRD(^,i[oiii]) is applied, i.e.. 



a(^,-^[Oiii]) 



/obs 

fobs X PRd(2, i[OIII]) 



(8) 



Here, /obs = 10,480/25,307 = 0.414 accounts for 
the incompleteness in the spectroscopic observa- 
tions of these "serendipity" targets, and is sim- 
ply the fraction of targets with observed spectra. 
The factor prd(-2, i[oiii]) accounts for the contri- 
bution of radio-weak type 2 quasars with optical 
SEDs similar to the radio-selected sources (and 
therefore would not be selected by the other algo- 
rithms), but which are too radio- faint to have been 
detected by the FIRST survey. We define it to be 
the probability that an object with a given redshift 
and [O III]5007 lumi nositv has a FIRST catalog 
match ([Becker et al ] [l995 'l within 2" of its optical 
position. In Wd.'Sl we show the type 2 quasar LF 
calculated both with and without this correction. 
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The distribution of radio luminosities in AGN 
and the relation, if any, between optical and ra- 
dio properties are the subjects of much debate 
(jivezic et al.l |2004| and references therein). Here, 
we determine Prd(z, ^[oiii]) from our data alone. 
We use the radio properties of type 2 quasars 
in our sample that are selected by at least one 
non-radio method (539 objects) to predict how 
many objects the Serendipity FIRST algorithm 
has missed because they fall below the flux limit of 
the FIRST survey. In the following calculation, we 
assume that the distribution of radio luminosities 
of a population of type 2 quasars is independent 
of their optical colors and morphologies, but may 
depend on their [OIII]5007 luminosities and red- 
shifts. 

We aim to find a simple analytic functional form 
for Pkd{z, i[oiii]) that best fits the available data. 
We define a statistic analogous to to assess the 
goodness of the fit. For a bin j in the space of 
redshift and [OIII]5007 luminosity, which contains 
A'^ objects with redshifts and luminosities Zi and 
Li {i = l...N) of which n are detected by FIRST, 
this statistic is given by 

where / = n/N and p = J^i PRD i^i, Li)/N. In the 
limit of large N and small bins (so that prd — con- 
stant within each bin), we show in the Appendix 
that u is distributed as with one degree of free- 
dom. Therefore, the central limit theorem can be 
applied to the sum of the values of u over all bins, 
and this value gives the probability that the given 
form of pp{d(-2, i[oiii]) tits the data. Importantly, 
this statistic is independent of the distribution of 
objects in the z-Ljom] plane; this distribution is 
strongly affected by selection effects. 

We expect that the probability of radio detec- 
tion decreases with redshift, since more distant ob- 
jects are dimmer, and increases with luminosity. 
We tried several functional forms for prd with this 
in mind and varied their parameters to minimize 
the sum of u over all bins. Our best-fit function is 
given by 

1 



1 



z/(0.15- 



O.l(logVi0-8.O)2)]2- 
(10) 

The comparison of observed FIRST detection 
rates with those calculated using this function 



is shown in Fig. |4l By visually examining the 
agreement between the observed detection rates 
and pRD (using slices and projections in z-Ljoiii] 
plane) we estimate that prd is determined to bet- 
ter than ±10% over most of the parameter space. 
Since Prd(-z, -^[oiii]) ^ 1 at low redshifts, as well 
as at high [OIII]5007 luminosities, the correction 
to the LF from this factor is not important in 
these regimes. In §5.31 we test whether the radio- 
detection rates of type 2 quasars derived in this 
section are statistically consistent with those of 
type 1 quasars. 

3.3. Results 

We derive the type 2 quasar luminosity func- 
tion for three ranges in redshift: z < 0.30, 0.30 < 
z < 0.50 and 0.50 < z < 0.83, using Eqs. [5]and 
m summed over 420, 175 and 145 objects, respec- 
tively. Results with and without correction for the 
probability of detection of radio-selected objects 
fi 33.2.4j) are shown in Fig. El 

The flattening of the 0.30 < z < 0.50 and 
0.50 < z < 0.83 LFs at low [OIII]5007 luminosities 
is due to the decreased sensitivity of all the tar- 
get algorithms to low-luminosity objects at high 
redshifts, mainly because of the algorithms' fiux 
limits (see also §3.4|) . This incompleteness is par- 
tially mitigated by the radio-detection correction, 
but even this result is still strictly a lower limit. 
Therefore, for these redshift ranges, the most use- 
ful information from the derived LFs is provided 
by the highest luminosity bins. The trend is re- 
versed for the z < 0.30 LF, where the most use- 
ful information is provided by the low luminosity 
bins, since these represent many detected objects 
(see also §0)) . 

While the high luminosity end of the LFs seem 
to suggest a tentative trend toward increasing 
number densities with increasing redshifts, we can- 
not make any conclusive statements because these 
are all lower limits. We return to the question of 
redshift evolution in the sample in ^5.2\ but it is 
difficult to disentangle the effects of incomplete- 
ness of the sample with that of redshift evolution, 
and we make no attempt to do so. 

Our result is consi stent with t he typ e 2 AGN 
luminosity function in lHao et al.l (|2005bl ). derived 
from objects with redshifts z < 0.30, selected from 
the SDSS Main Galaxy spectroscopic database. 
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Figure |6] shows the high-luminosity end of this LF. 
We find good agreement with our derived type 2 
quasar LF in the redshift and luminosity ranges 
in which the two sa mples over l ap. In this Figure, 
the luminosities of iHao et al. I (l2005lj) have been 
shifted up by 0.14 dex to account for the difi'er- 
ence in their spectrophotometric flux cahbration 
(see §2.3p . A more r ecent measurement o f the type 
1 AGN Ha LF by ICreene fc Hoi (|2007h is ~ 40 
ti mes lower a t low lu minosities than that derived 
bv lHao et al. I (l2005bl ) because of their more strin- 
gent sample selection criteria, but the two LFs are 
consistent with each other at the high-luminosity 
end (iv[oiii] '■'^ IO^Lq), that is, at luminosities rel- 
evant to our comparison. 

3.4. How Many Type 2 Quasars have we 
Missed? 

In principle, if we knew the distribution of op- 
tical SEDs and radio properties of type 2 quasars, 
we could calculate the probability with which each 
algorithm would select a type 2 quasar of a given 
[OIII]5007 luminosity at every redshift. We could 
then use this information to correct the luminosity 
function for objects that are not in our sample be- 
cause they do not have the optical colors, apparent 
magnitudes, optical morphologies, or radio fluxes 
that are targeted by the SDSS algorithms. The 
optical spectrum of a type 2 quasar is the sum of 
the host galaxy spectrum, the n arrow lines and 
scattered light from the AGN (IZakamska et al 



200,4 l200fih . None of these components can be ne 
glected, the relative strengths of the components 
vary from object to object, and even the shape 
of the host galaxy continuum matters for selec- 
tion. For example, the Balmer break of the host 
galaxy moving across the r and i fllters at red- 
shift ^ 0.35 produces a difference in colors suffi- 
cient to enable or disable selection by the Low- 
z QSO algorithm in several objects. This se- 
lection effect results in the weak feature in the 
[OIII]5007 luminosity-redshift distribution seen in 
upper right-hand panel of Fig. |3l 

With the above complications in mind, we 
do not attempt to model the spectra of type 2 
quasars. Instead, we derive quantitative estimates 
of the contributory power of the various target 
algorithms to the selection of type 2 quasars, di- 
rectly from the data. We do not use these results 
to apply corrections to the LF, in keeping with our 



'lower limit' approach; but they are nonetheless 
useful in guiding the interpretation of the derived 
luminosity functions. 

We define the contributory power of a given tar- 
get algorithm to be the fraction of objects in the 
sample (871 objects with £[oiii]> lO^'^L©) that 
we would select if we use this algorithm alone. We 
determine what this fraction would be if all objects 
were placed at a certain redshift (following ij3.2p . 
Figure [7] shows the resulting function for assumed 
redshifts in the range of 0.1 to 0.8 (in steps of 
0.01) and for three ranges in [OIII]5007 luminos- 
ity. Not surprisingly, we find that the algorithms 
with deeper magnitude limits (Serendipity FIRST 
and High-z QSO) do better at selecting objects at 
high redshifts than do the shallow ones (Galaxy 
and Low-z QSO). The color-based selection crite- 
ria of the Low-z QSO and High-z QSO algorithms 
cause the peaks and dips in their curves. We find 
that all four algorithms are poor at selecting low 
L[Oiii] objects at high redshifts, which is why the 
derived type 2 quasar LF shows evidence for in- 
completeness in this region. 

3.4-1- Contributory Power of Color-Based Tar- 
get Algorithms 

An alternative approach to estimate the con- 
tributory power of the color-based target algo- 
rithms to the selection of type 2 quasars is to ex- 
amine their overlap with the non-color based algo- 
rithms. Of the 258 Galaxy targets with redshifts 
z < 0.2, 254 (98%) are brighter than the magni- 
tude limit of the Low-z QSO algorithm (jpsF = 
19.1) but only 12 of these (5%) have colors that 
satisfy the Low-z QSO selection criteria. Between 
z = 0.2 and 0.4, there are 137 Galaxy targets; 113 
of them (82%) have ipsF < 19.1, but only 1 object 
(0.9%) was selected by the Low-z QSO algorithm. 
The small amount of overlap between the Low-z 
QSO and Galaxy targets indicates that the Low- 
z QSO algorithm is missing many type 2 quasars 
because it is only sensitive to a limited region of 
color space. There is no overlap at all between the 
High-z QSO and Galaxy algorithms. The colors of 
many type 2 quasars are dominated by the colors 
of their host galaxies, especially for lower [OIII] 
luminosities, and neither Low-z QSO nor High-z 
QSO algorithms target such objects. 

Out of the 275 objects selected by the Serendip- 
ity FIRST target algorithm, 74 have ipsp < 19.1 
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and 15 of these (20%) are targeted by the Low-z 
QSO algorithm; 245 have ipsF < 20.2 and only 
12 of these (5%) are targeted by the High-z QSO 
algorithm. These comparisons suggest that over 
the full redshift range, the Low-z QSO and High-z 
QSO algorithms alone select only 20% or less of 
type 2 quasars. 

3.4-2. Low-redshift Objects and the Galaxy Tar- 
get Algorithm 

In this section, we discuss the interpretation of 
our LF calculation in the lowest redshift range, 
z < 0.30. We see from Fig. [3] that there are no 
objects in our sample with i[oiii]> lO^'^L© found 
at these redshifts, and only 1 out of 17 objects with 
-^[oiii]> W^'^Lq is at z < 0.5. At face value, the 
scarcity of detected objects may be interpreted as 
a significant drop in the ratio of number densities 
of type 2 quasars at these luminosities. However, 
we argue here that this may instead be due to a 
combination of the selection effects and the small 
volume covered by the lowest redshift bin. 

Most type 2 quasars in the redshift range z < 
0.3 were selected by the Galaxy algorithm. The ef- 
ficiency of the Galaxy algorithm begins to drop at 
redshifts around z = 0.2 due to its shallow mag- 
nitude limit (Figs. [3] and [7|). Quantitatively, an 
object with L[oiii]= IQ^'^Lq at z = 0.3 with a 
rest-frame equivalent width of 500A (correspond- 
ing to the median of the i[oiii]/EW[oiii] distribu- 
tion for our sample; Paper I) would have an r band 
AB magnitude of around 17.8 mag, and would 
therefore not be selected by the Galaxy algorithm 
(which has a limiting magnitude of 17.77 mag). 
As we discussed in the previous section, these ob- 
jects may be missed by the color-based target al- 
gorithms because they do not lie in the specific 
regions of color-color space these algorithms were 
designed to target. The radio-based Serendipity 
FIRST algorithm would have a chance of iden- 
tifying high-luminosity, low-redshift objects that 
happen to be radio-bright, but even the luminous 
objects at low redshift may fall below the radio 
fiux limits of the FIRST survey. More impor- 
tantly, only - 40% of Serendipity FIRST targets 
have observed spectra ( ^3.2.4p . Therefore, the fact 
that we have no type 2 quasars with z < 0.30 and 
i[Oiii] > 1O^'^L0 may be entirely due to the incom- 
pleteness of SDSS spectroscopic target selection. 



4. Type 1 Quasar Luminosity Function 

In order to put our type 2 quasar LF into con- 
text, we need to compare it with the type 1 quasar 
LF at similar redshifts and luminosities. In this 
section, we derive the l/Vma^ [OIII]5007 LF from 
a complete sample of type 1 quasars using the 
formalism we have applied to the type 2 quasar 
sample f H3.ip . We discuss the sample selection 
in H4.ll measurement of [OIII]5007 luminosities in 
^4^^ the calculation of the LF and results in i )4.3[ 
and comparison with previous work in §4.41 The 
type 1 LF for the three redshift ranges is shown 
in Fig. [11] In the following section, ^jSJ we com- 
bine this result with the derived type 2 quasar LF 
to calculate the ratio of space densities of the two 
populations. 

4.1. Type 1 Quasar Sample 

We select a complete sample of 8,003 type 
1 quasars from the S PSS DR5 Quasar Catalog 



(jSchneider et al.M2007l ). The catalog consists of 



77,429 spectroscopically-confirmed quasars in the 
redshift range 0.08-5.41, with absolute magni- 
tude in the i band. Mi < —22 mag and full 
width at half-maximum (FWHM) of lines from 
the broad-line region greater than 1000 km s^^. 
From this catalog, 31,999 objects are targeted 
with the final version of the Low-z QSO algorithm 
(iRichards et al...2 002'). Of these, 8,003 objects sat- 
isfy the redshift cut z < 0.83 and are included in 
our sample. The effective area of our sample is 
4,041 deg^, calculated using the same procedure 
as in 

The redshift and [OIII]5007 luminosities of 
these objects are shown in the bottom-right panel 
of Fig. [3l There is no explicit lower cut in red- 
shift, but there are very few objects with redshifts 
z < 0.15 because of the strong redshift evolution 
in the type 1 quasar population and the bright 
limit of the Low-z QSO target algorithm. 

4.2. Measurement of [QUI] 5007 Luminosi- 
ties 

In the calculation of the type 1 quasar LF, we 
use [OIII]5007 luminosities measured from Gaus- 
sian fits, after careful subtraction of Fell emission. 
The average contribution of Fell emission to the 
total flux over the wavelength range 5007A± 1.7x 
FWHM is ^ 9%. We exclude the region con- 
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taining H/3 and [OIII]4959,5007 and find the best 
fit to the spectrum in the form of a p ower law 
plus the Fell template o fH oroson fc Gre en (1992). 
We then fit a set of four Gaussians to [OIII]4959, 
[OIII]5007, and narrow and broad components of 
H/3. 5,297 quasars (67% of the sample) have 
-^[Oiii]> lO^-^i©, and 221 quasars (2.8% of the 
sample) have L[oiii]> lO^-^L©. 

We also obtain a non-parametric line lumi- 
nosity by integrating the detected flux density 
and subtracting the contribution from the con- 
tinuum, Fell, H/3, and [OIII]4959. For most ob- 
jects, the Gaussian and integrated line fluxes are 
in good agreement. Outliers from the integrated 
vs. Gaussian locus were visually inspected, and 
187 spectra (2.3% of the sample) for which the fit 
failed (due to low signal-to-noise ratio, bad pixels, 
etc.) were discarded. 

4.3. l/VJnax Luminosity Function 

Objects in the type 1 quasar sample were all 
selected by the Low-z QSO target algorithm, de- 
scribed in i i3.2.2l In the case of type 2 quasars, we 
considered both the magnitude limits and color- 
based selection criteria. In contrast, type 1 quasar 
colors d o not evolve strong l y wit h redshift for 
z < 2.2 ([Richards et al.ll200ll [200^ . so we do not 
expect the color-based selection criteria to be im- 
portant in this case. For each object, we deter- 
mine the available comoving volume Knax based 
only on the faint magnitude limit of the survey, 
ipSF = 19.1, since the volume corresponding to 
the bright limit of the survey is negligibly small. 

As in H3.2[ we calculate scaled and k-corrected 
PSF magnitudes as a function of assumed redshift 
from the observed spectrum of each object. We 
determine ^max, the redshift at which the PSF 
magnitude reaches the limit of the survey. The 
available volume Vmax is then given by 



V 
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where fl = 4041 deg^ is the effective survey area 
for this sample, Zi and Z2 are the edges of the 
redshift range for which we are calculating the LF 
and {dVc/dz)dz is the comoving volume element 
in the redshift interval dz. 

Figure [S] shows the 1/Knax [OIII]5007 luminos- 
ity function for type 1 quasars (calculated using 



Eqs. [5] and |6]), for the same redshift ranges as 
for type 2 quasars. Our results reflect the posi- 
tive red shift evolution of the type 1 quasar pop- 
ulation ( Richards et al. 2006[) . The turn-over of 



the 0.50 < z < 0.83 LF at low luminosities is an 
artifact due to the difficulty of measuring weak 
[OIII]5007 lines at these redshifts, since the line 
falls on the red end of the observed spectrum 
where the signal tends to be noisier. 

4.4. Comparison with the Broadband Type 
1 Quasar LF 

In this section, we test whether our measured 
[OIII]5007 type 1 quasar LF is consistent with the 
observed broadband LF, based on type 1 quasars 
from the SDSS DR3 quasar catalog, with redshift s 
0.30 < z < 0.68 (Table 6 of lRichards et al.ll2006[) . 
The measure of luminosity used there is the con- 
tinuum luminosity around rest-frame wavelength 
25OO A, corresponding to t he SDSS i band at z = 2 
(e.g., iBlanton et al. 2003 ). expressed as a broad- 
band absolute magnitude M25oo- This test serves 
two purposes. First, it is a sanity check for our 
LF calculation. Second, it probes the correlation 
between [OIII]5007 line luminosity and continuum 
luminosity in type 1 AGN. This is important for 
justifying our assumption that [OIII]5007 luminos- 
ity traces bolometric luminosity sufficiently well 
for our purposes (see ^6.ip . As expected, we find 
a strong correlation between M2500 and [OIII]5007 
luminosity in our type 1 quasar sample (Fig. [9]) . 
The mean relation is 

-^[OIII] \ 

Lq ) 

consistent with a linear relation, with a slope of 
0.95 in the logL25oo^ogL[oiii] plane. The scatter 
in L[oiii] at fixed Af25oo is consistent with a log- 
normal distribution of width 0.36 dex. 

To convert the broadband LF, <&(M25oo), into 
an [OIII]5007 LF, $conv(i[oiii])j convolve the 
broadband LF with the mean i[oiii]--^^2500 rela- 
tion Eq. [T2I with a log-normal scatter of width 
a = 0.36 dex, i.e., 



log 



-0.38M; 



2500 



0.62, 



(12) 



*co„v(/:) - / $(M)exp 



(£ - C{M)f 



2a^ 



dM, 
(13) 

where £ = log(L[oiii]/io) and M = Ahsoo- Fig- 
ure [TO] shows that there is a good agreement be- 
tween the converted LF (short-dashed curve) , and 
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our measured [OIII]5007 LF for 0.30 < z < 0.50 
(solid line, open circles). The turn-over at the low 
i[Oiii] end of the converted LF is due to the ab- 
solute magnitude cut-ofF, since there are no ob- 
jects fainter than M2500 ^ —23 mag to scatter 
into these ipiii] bins. 

The shape of the converted LF is sensitive to 
the assumed slope and scatter in the i[oiii]~^'jf2500 
relation. In particular, if zero scatter in the i[oiii]- 
relation is assumed, the resulting converted 
[OIII]5007 LF is inconsistent with our derived LF 
(long-dashed curve). With scatter properly taken 
into account, we find that the two independent 
derivations of the LF match. This suggests that 
we can successfully convert an [OIII]5007 LF to a 
broadband LF (and vice versa). 

5. Comparison of Type 1 and Type 2 
QuasEir Luminosity Functions 

In this section, we combine the results of the 
previous sections to compare the type 1 and type 
2 quasar populations. Most importantly, we can 
now use the derived [OIII]5007 type 1 and type 2 
quasar luminosity functions to constrain the frac- 
tion of type 2 quasars as a function of [OIII]5007 
luminosity. We discuss the calculation, results and 
their interpretation in ^15. II We compare the two 
populations in terms of their redshift evolution in 
and their radio properties in 



5.1. Type 2 Quasar Fraction 

Figure [TT] shows that the [OIII]5007 type 1 and 
type 2 quasar luminosity functions are comparable 
in the regimes where we expect our type 2 quasar 
LF to be least affected by selection bias, i.e., the 
low-luminosity regime (< lO^L©) for the redshift 
range z < 0.30 and the high-luminosity regime 
for redshift ranges 0.30 < z < 0.50 and 0.50 < 
z < 0.83. From these two functions, we directly 
calculate the type 2 quasar fraction, i.e., the ratio 
of type 2 to total (type 1 -I- type 2) quasar number 
densities. We determine the number density of 
quasars in a given luminosity range by integrating 
the LF over that range. For this calculation, we 
use our best lower bound to the type 2 quasar LF 
(which includes a correction for the probability of 
radio detection, i j3.2.4p . 

Figure [12] shows the calculated type 2 quasar 
fraction as a function of [OIII]5007 luminosity. For 



the z < 0.30 bin, we find that the type 2 quasar 
fraction is about 60% for L[oiii]= lO* '^ to lO^i©. 
For z < 0.30 and 0.30 < z < 0.50, our strongest 
limits come from the highest luminosity bin, which 
constrain the type 2 quasar fraction to be at least 
40% and 60% at i[oiii] ~ 5 x lO^L©, respectively. 

Our derived type 2 quasar fraction is a lower 
limit for the following reasons: (i) our derived type 
2 LF is a lower limit (see beginning of @; (ii) the 
combination of selection effects and redshift evo- 
lution artificially lowers the type 2 quasar fraction 
that we derive (as we describe in §5.2^ : (iii) there 
are indications that the [OIII]5007 line is slightly 
more extinctcd in type 2 quasars than in type 1 
quasars (t j5.3l and and (iv) we have assumed 

that the [OIII]5007 luminosity is independent of 
the obscured fraction and serves as a tracer of 
the bolometric luminosity. In practice, it is likely 
that the [OIIIJ5007 luminosity is higher for objects 
with a larger opening angle, since in this case more 
material is illumi nated by the central source. As 
demonstrated by KrolikI ( 1999t ). this dependence 
leads to a bias favoring unobscured objects in any 
[OIII]-flux selected sample. 

There has been a substantial amount of work to 
determine the obscured quasar fraction at differ- 
ent luminosity and redshift regimes. Some of these 
results are shown together with ours in Fig. [121 
Nearest to our approach are previous determina- 
tions from emission-line luminosity functions of 
low-luminosity (type 1 and t ype 2) AGN in the 
SDSS, at redshifts z - 0.1. iHao et al.l (flooHE) 
found that type 2 AGN make up about 60% of the 
AGN population at I/[nTTT] ^ lO^L© , and about 
30% at ~ 3 X lO'^L©. ISimpsonI (|2005') also clearly 
finds the decreasing trend in type 2 fraction with 
luminosity, but found substantially higher type 2 
fractions (open diamonds). 

Determinations from hard X-ray (2-10 keV) se- 
lected samples, at redshifts z ^ 3, suggest that 
the fraction of obscured quasars is large at low 
luminositie s, and then decreases at higher lu- 
minosities (lUeda et al.ll2003l: ISzokolv et al.ll2004 



Bargcr et al.ll2005l:lMar kwardt ct al. 2005; Trcis ter et al 
i2006 : .Beckmann et al.L2006i : .Sazonov et al ,,200j 
Figure [m shows results from Ueda et. al. (2003; 
open squares) and Hasinger (2008; open circles 
and dashed line). We have converted the hard 
(2-10 keV) X-ray luminosities to [OIII]5007 lu- 
minosities by shifting by 1.59 dex (toward lower 
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luminosities) following iHeckman et alj (j2005l) . X- 
ray surveys may be missing Compton-thick ob- 
jects which might constitute a significan t fra c- 
tion of ah AGN (iHeckmanet alJ [20051: IPolletta 
2006t Martfnez-Sansigre et al. 2007 ). so these val- 
ues should also be treated as lower limits to the 
true obscured quasar fraction. Also shown are re- 
sults from radia data (Grimes et al. 2004; open 
triangles) . 

Fro m samples selected from infrar ed observa- 



tions, Martmez-Sansigre et al 



(|2006l) estimated 
the fraction of obscured quasars at z ~ 2 to be 
~70%, with a large uncerta i nty b ecause of small 



sample size. iPolletta et al.l ()2008[ ) estimated the 
fraction of obscured quasars to be 60 — 65% at 
bolometric luminosities 10'*^~'*^erg s~^ {L[om]^ 
IO^Lq), using a sample of type 2 quasars with 
Si absorption. Again, these values should be 
taken as lower limits because, e.g., the former 
sample excludes radio- weak objects, while the lat- 
ter excludes sourc es with featureless IR spectra 
(|Sturm et al.ll2006l ). 



Recently, iTreister et al.l (|2008[ ) determined the 
type 2 quasar fraction by calculating the fraction 
of the total light emitted by obscuring material in 
the IR in type 1 AGN (filled triangles). Shown are 
the results assuming the parameter /12 — 0.06, the 
fraction of the total dust-reprocessed luminosity 
falling within the MIPS band (11-13 in the rest 
frame at z 1), and with model-dependent error 
bars. We have converted their quoted bolometric 
luminosities to [QUI] 5007 luminosities using the 
approximate conversion in §6.11 

Overall, our lower limits on the type 2 quasar 
fraction are consistent with previous determina- 
tions. In practice, however, we suspect that the 
true type 2 quasar fraction is significantly higher 
than our lower limits, for the reasons discussed 
above and in Sj6l This simply underscores the 
need to treat most determinations of the obscured 
quasar fraction as lower limits, since a substan- 
tial part of the population may be missed by 
wavelength- and method-specific selection criteria. 

5.2. Comparison of Redshift Evolution 

In this sec tion, we apply the V/Vmux test 



(|Schmidtill968[ ) to probe redshift evolution within 
the type 1 and type 2 samples. For a uniformly- 
selected, non-evolving population of objects. 



V/Vmax is expected to be distributed uniformly 
between and 1 and to have a mean value of 
0.5. For an object with redshift within the range 

(zi,Z2), V/Vn^ax IS givcu by 



(14) 



where Zobs is the observed redshift of the object 
and (dVc/dz)dz is the comoving volume element 
in the redshift interval dz. If there is positive red- 
shift evolution within the redshift range, the mean 
values of T^/Vmax over the sample are greater than 
0.5. Incompleteness in the sample can either mask 
or enhance this effect, by weighting the sample to- 
ward the lower or higher end of the redshift range. 

In Table [31 we list the mean values of F/Vmax 
calculated from the type 1 and type 2 quasar sam- 
ples for the three redshift ranges we used pre- 
viously. Our type 1 quasar sample is complete, 
and the (F/Vmax) values greater than 0.5 indicate 
the positive redshift e volution within the sample 
( Richards et al.ir2006 ). On the other hand, our 
type 2 quasar sample is incomplete, so the devi- 
ation of (V^/Knax) from 0.5 is due to a combina- 
tion of redshift evolution and selection effects (see 
' i'6A\ . which we do not attempt to disentangle. For 
two out of the three redshift ranges, (l^/Vmax) val- 
ues are lower for the type 2 quasar sample than for 
the corresponding type 1 quasar sample. 

Since our derived LFs are calculated for broad 
ranges in redshift, they are a weighted average of 
the true LFs over the redshifts within this range. 
The above results suggest that the type 1 quasar 
LFs are weighted toward the higher end of the 
redshift range, and therefore toward higher num- 
ber densities, relative to the type 2 quasar LFs. 
Therefore, the combination of selection effects and 
redshift evolution of quasar number density arti- 
ficially lowers the type 2 to type 1 ratio that we 
derive from these LFs. 

5.3. Comparison of Radio Properties 

In this section, we test whether the radio prop- 
erties of type 2 quasars are statistically consistent 
with those of type 1 quasars. We test this hy- 
pothesis to the same limited extent that we have 
employed in t ^3.2.4| - that is, we test whether non- 
radio-selected type 2 and type 1 quasars have the 
same detection rate with FIRST within a match- 



14 



ing radius of 2". We use a variation of the statistic 
defined in Eq. |9] for comparing two populations of 
objects: 

{fi-f2? N,N2 
^h{l-h)f2{l-h)Ni+N2 

Here A^i and N2 are the number of type 1 and type 
2 quasars, respectively, in a given bin in the z- 
L[Oiii] plane and /i and /2 are the fractions of ob- 
jects in this bin that are detected by FIRST. The 
distribution of the statistic described by_^£;_^ 
(insp ired by two-sample statistics from IPeacock 
19831 ) is not formally independent of fi and Ni in 
the same sense that the distribution of Eq.[9]is in- 
dependent of A'' and p (as shown in the Appendix). 
Using Monte-Carlo simulations we found that the 
distribution of the value in Eq. [HI is very close to 
that of with one degree of freedom, as long as 
/i and /2 are not too close to or 1 and as long 
as the number of objects in each bin is sufficiently 
large {Ni > 10). 

We find that the radio-detection rates of the 
two populations are statistically indistinguishable 
for all except the low-redshift, low-luminosity ob- 
jects. Type 2 quasars with z < 0.25, 8.0 < 
log(L[oiii]/-^0) < 8.5 are significantly more likely 
to be detected in FIRST than type 1 quasars. One 
possibility is that the [OIII]5007 emission line is 
systematically more extincted by dust in type 2 
quasars than in type 1 quasars. Correcting for 
0.1 dex of this putative extinction removes the de- 
tection rate difference everywhere in the z-L[oiii] 
plane. 

6. Discussion 

We have presented our main results in terms 
of [OIII]5007 luminosity, a proxy for AGN activ- 
ity that is directly measurable from the observed 
spectra. In this section, we discuss two issues that 
are important for the interpretation of these re- 
sults: the relation between [OIII]5007 and bolo- 
metric luminosity and the amount of extinction of 
the [OIII]5007 fine (gUI] and 1121 respectively). We 
conclude that despite these effects, our approach 
to calculate lower limits on the type 2 quasar lu- 
minosity function and the type 2 quasar fraction 
remains valid. In ii6.3[ we discuss implications of 
our results for estimates of the accretion efficiency 
of supermassive black holes using measurements 
of AGN space densities. 



6.1. [OIII]5007 Luminosity as Tracer of 
Bolometric Luminosity 

The [OIII]5007 emission line luminosity is ar- 
guably the best available proxy for AGN activ- 
ity for optically-selected obscured quasars. This 
line is emitted by the narrow-line region, which 
extends outside the obscuring material thought 
to surround the broad-line region of typ e 2 AGN 
(|Antonuccilll993t lUrrv fc Padovanilll995l) . Obser- 
vational support for this assumption comes from 
the similarity of I R/[OIII] ratios in type 1 and type 
2 AGN found bv lMulchaevI (|l994[ ). Moreover, the 
[OIII]5007 line has the advantage of beii ig strong 
and e asily detected in most AGN, and ISimpson 
( 19981 ) and iKauffmann et al" ( 2003 ) have shown 
that it is a good tracer of AGN activity and is not 
severely contaminated by star formation. 

For our type 1 quasar sample, we found that 
[OIII]5007 luminosity correlates strongly with 
broadband luminosity at around 2500 A (Fig. [9]), 
as well as with the monochromatic continuum lu- 
minosity measured at rest-frame 5100A (linear 
correlation coefficient = 0.53). More importantly, 
we have demonstrated in ii4.4l that given the ob- 
served mean and scatter of the correlation between 
[OIII]5007 and continuum luminosity allows to 
quantitatively convert the type 1 quasar broad- 
band LF to a [OIII]5007 LF that is consistent 
with our independent measurement. 

Of course, [OIII]5007 luminosity is not a per- 
fect tracer of bolometric luminosity. Indeed, there 
is substantial scatter between [OIII]5007 and con- 
tinuum luminosity for type 1 quasars (see Fig. (9] 
Netzer et al. 20061 and references therein). It is 
plausible that this scatter reflects a real, physical 
difference in covering fractions of the narrow-line 
region. As mentioned in ^5.1\ this would imply 
that our [OIII]-flux selected sample is biased to- 
ward type 1 sources and would artificially lower 
our derived type 2 quasar fraction. This is in line 
with our approach to calculate lower limits on this 
quantity. 

Rough estimates of bolometric luminosities can 
be found by, first, using the Af25oo — -^[oiii] relation 
(Eq.[T2), and then the well-measured average type 
1 quasar spectral energy distributions. For exam- 
ple, an [OIII]5007 luminosity of 3 x lO^L© corre- 
sponds to an intrinsic UV luminosity of M2500 ~ 
—26.6 mag. Next, we use the average quasar SED 
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from IVanden Berk et alj (|200lh to calculate the 
corresponding luminosity in the B ba nd. Finally; 



we use bolometric corrections from Marconi et al 



(|2004l ). derived from a template AGN spectrum 
and with model 1-sigma uncertainties of ~ 0.05 
dex. We estimate that the corresponding bolomet- 
ric luminosity is ^ 1 x lO**^ erg s^^ or 3 x lO^'^i©. 

6.2. Reddening and Extinction 

The [OIII]5007 line, which we use as a proxy for 
AGN activity, is emitted by an extended narrow- 
line region. It is expected to be affected by extinc- 
tion due to interstellar dust, located either within 
the narrow-line region itself or in the intervening 
interstellar matter of the host galaxy. In this sec- 
tion, we estimate the magnitude of this extinction 
for type 2 quasars and discuss its effects on our 
measurements of the luminosity functions and the 
type 2 quasar fraction. 

We attempt to estimate narrow-line region ex- 
tinction for the luminous subsample of type 2 
quasars (with L[oiii]> lO^i©) by determining 
their Balmer line ratios. We measure Ha/H/3 for 
44 objects (with z < 0.4, for which the Ha line is 
observed) and II/3/II7 for 203 objects (for which 
the H7 line has sufficiently high S/N). We use 
fluxes measured using the non-parametric line fit- 
ting procedure described near the end of §2.31 and 
illustrated in Fig. [1] 

Our results are shown in Fig. 1131 The left panel 
shows II/3/II7 vs. [OIII]5007 luminosity. The me- 
dian value of II/3/H7 is 3.06 (dotted line), which 
corresponds to an extinction of [OIII]5007 of 2.6 
mag for a Milky way extinction curve with i? = 3.1 
and a foreground obscuring screen. Error bars 
are based on the difference between the flux mea- 
sured from our non-parametric fitting procedure 
and that derived from a Gaussian fitting procedure 
or from integration over the emission line (after 
subtraction of the continuum). The right panel 
shows Ha/H/3 vs. H/3/H7 for the subset of ob- 
jects with both measurements. The arrow shows 
the extinction vector for the Milky Way extinction 
curve, whose length corresponds to an [OIII]5007 
extinction of 2.2 mag. Other commonly used red- 
dening laws (e.g., the Small Magellanic Cloud ex- 
tinction curve) yield very similar extinction vec- 
tors. The data are clearly inconsistent with stan- 
dard reddening laws. For example, the median 
Ha/H/3 = 4.25 corresponds to an extinction of 



1.2 mag, which is '-^ 1.4 mag lower than that de- 
rived from the median H/3/H7. There is no sta- 
tistical evidence for a correlation between Ha/H/3 
and H/3/H7 (correlation coefficient = 0.3). Fur- 
thermore, the data are inconsistent with a linear 
regression from the case A-B values, and there- 
fore cannot be described by a dust screen model, 
no matter what the extinction law is. 

Early studies of Balmer line intensities of 
quasars have indicated that these cannot be ex- 



plained by standard recombination theory and 
a standard dust reddenine law (1 Anderson 1970t 


Adams et al.lll975: Ba 


dwinll975HOsterbrock et al. 


19761: Osterbrockll977 


). Other processes that may 



play an important role include collisional excita- 
tion from higher energy levels, self-absorption of 
the Balmer lines, continuum optical depth, scat- 
tering and fluorescence ( tC apriotti, 1964a . b: Netzeil 



19751 iKrolik fc McKeel[l978l ). Given the complica 



tions, we do not apply dust extinction corrections 
to our [OIII]5007 luminosity measurements. 

Our use of extinction-uncorrected line luminosi- 
ties means that in our calculation of the luminosity 
function, objects with [OIII]5007 lines affected by 
extinction have been assigned to fainter luminos- 
ity bins. Hence, the LF has been shifted toward 
fainter luminosities, or equivalently, toward lower 
space densities, relative to an extinction-corrected 
LF. Moreover, if line extinction is stronger in type 
2 than in type 1 sources, then the type 2 quasar 
fraction would also be underestimated. Several in- 
dependent measure ments have suggest ed that this 
is indeed the case ( Netzer et al. 20061 and refer- 



ences there in) . Reddening estimates from Balmer 
decrements suggest that the narrow-line regions 
of type 2 sources are m ore heavily reddened than 
those of type 1 sources ( De Zotti fc Gaskell 19851 : 



Dahari &: De Robertid 19881 ). Recent IR observa- 
tions strongly indicate that the narrow-line regions 
of type 2 quasars are significantly more obscured 
than those of type 1 q u asars, by up to a factor 
of 10 (|Haas et al.ll2005l ). iNetzer et all ([2006 ) sug- 
gest that [OIIIJ5007 luminosity is a factor of '--^ 2 
larger in type 1 than in type 2 quasars for the same 
2-10 keV X-ray luminosity, as would be the case 
if type 2 sources are more reddened than type 1 
sources. Finally, such a trend is also hinted at by 
our analysis of radio properties of type 1 and type 
2 sources ( §5.3^ . Therefore, the use of luminosi- 
ties uncorrected for extinction is consistent with 
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our approach to calculate lower limits on both the 
type 2 quasar LF and the type 2 quasar fraction. 

6.3. Implications for the Black Hole Mass 
Function and Accretion Efficiency 

If black holes grow primarily by accretion, the 
total luminosity emitted by accretion processes 
over the lifetime of the Universe is directly re- 
lated to the accum ulated mass of local super- 
mass ive black holes ( Soltan 1982 : Yu fc Tremaind 
Recent measurements of both the mass 
density of supermassive black holes and the to- 
tal luminosity emitted by matter accreted onto 
them can be reconciled if the radiative efficiency 
(the fraction of the accreted mass that is con- 
verted directly into obser ved quasar radiation) 
is quite high e > 0.1 (IYu fc Tremaind 12002 : 



Marconi et all I2OO4I: [Barger et alJ l2005h . In this 



section we discuss qualitatively how these calcula- 
tions are affected by the existence of a large pop- 
ulation of obscured quasars. 

What matters for quantifying the growth of the 
black hole is its total energetic output ('bolomet- 
ric luminosity' Lboi), since that is the value that 
can be related to the accreted mass through radia- 
tive efficiency. Some of the optical, UV and X-ray 
radiation is intercepted by the obscuring material 
and then re-emitted in the IR more or less isotrop- 
ically, but no new energy is generated. The frac- 
tion of energy reprocessed this way roughly equals 
f{Lboi)Lboh where /(Lboi) is the fraction of ob- 
scured AGN in the population and ATTf{Lboi) is 
the solid angle covered by obscuring material as 
seen from the central engine. Since the IR emis- 
sion is directly related to the obscuration frac- 
tion, it is possible to derive the obscured fraction 
from the IR-to-optic al ratio in unobscured quasars 
(|Treister et al.ll2008l ). 

Some of the accretion efficiency calculations are 
based on the optical LF of type 1 quasars (e.g., 
Yu fc Tremaind I2OO2 ). In this case, the luminos- 



ity of each individual quasar is overestimated by 
including the IR emission in the bolometric cor- 
rection, but the total number of objects is under- 
estimated. The net result is that the contribu- 
tion of quasars with bolometric luminosities Lboi 
to the luminosity budget should be augmented 
by roughly 1/(1 — P{Lboi)) to account for the 
obscured sources. Other recent calculations are 
based on the hard X-ray luminosity functions of 



AGN (jMarconi et al.ll2004 iBarger et all 120051) in 

an effort to include both obscured and unobscured 
AGN, but this procedure misses the contribution 
from Compton-thick sources. The most recent 
accretion efficiency estimates are based on com- 
bining optical and X-ray da ta to produce bolo- 
metric luminosity fu nctions (jHopkins et al.l 12007 : 



Shankar et al.ll2007l ). Nevertheless, the fraction of 



Compton-thick sources which would be missed by 
both X-ray surveys and type 1 AGN surveys re- 
mains a major uncertainty in these methods. If 
fc{Lboi) is the Compton-thick fraction of type 2 
AGN, then the contribution of AGN with bolo- 
metric luminosities Lboi should be increased by 
1/(1 - fn(Lb oi)f(Lboi)). Assuming fc > 0.5 



(|Risalitilll999[ ) and /> 0.5 (this work) for quasars 
around and above the luminosity function break, 
we find that both types of methods (those based 
on the optically selected AGN and those based on 
X-ray selected AGN) may underestimate the AGN 
luminosity budget by 30% or more. This means 
that the calculated accretion efficiency would be 
underestimated by the same amount. 

7. Summary and Conclusions 

In this paper, we present the largest sample 
of type 2 quasars to date. Table [T] lists the 887 
optically-selected sources with redshifts z < 0.83. 
These objects are selected from the spectroscopic 
database of the SDSS on the basis of their emission 
line properties. Candidate sources are required to 
have no broad (FWIIM> 1100 km/sec) compo- 
nents in their permitted emission lines. To dis- 
tinguish type 2 AGN from star-forming galaxies, 
we use the standard line diagnostic diagrams in- 
volving [0III]/H/3, [Nil] /Ha and [SII]/Ha line ra- 
tios (at redshifts z < 0.3) or require other signs 
of a hidden AGN, such as the presence of the 
[NeV]3346,3426 emission hues (at redshifts 0.3 < 
z < 0.8). We place a lower limit on the [OIII]5007 
line luminosity of lO^'^L©, ensuring that the bolo- 
metric luminosities of these objects are above the 
classical Seyfert/quasar separation of 10^^ erg s~^. 

For this sample, we calculate the [OIII]5007 
luminosity function using the 1/Vinax method 
(Fig. [6]). The selection function that we calculate 
for each object involves taking into account four 
different spectroscopic target selection algorithms, 
of which two (Low-z QSO and High-z QSO) are 
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color-based, one (Galaxy) is based on optical mor- 
phology and one (Serendipity FIRST) is based on 
the radio properties. By comparing how well each 
algorithm performs at selecting type 2 quasars, 
we find that color-based algorithms are rather in- 
effective and select no more than 20% of objects. 
Indeed, there is no well-defined region of the color- 
color space where type 2 quasars concentrate, even 
when only narrow redshift ranges are considered. 
The success of the SDSS in selecting a large num- 
ber of type 2 quasars is due to the unprecedented 
size of the spectroscopic survey and the multitude 
of different target selection algorithms which allow 
for serendipitous objects. 

We extend the [OIII]5007 AGN luminosity 
function to luminosities about 2 orders of mag- 
nitude higher than was previously done, up to 
i[oni] — 10^'^ Lq. This value corresponds to an 
intrinsic UV luminosity of M2500 = — 28 mag and 
a bolometric luminosity of 4 x 10"*^ erg s^^. 

We also derive the [QUI] 5007 luminosity func- 
tion for a complete sample of 8,003 z < 0.83 type 
1 quasars taken from the SDSS quasar catalog and 
find it to be in excellent agreement with other mea- 
surements. We can then directly compare the lu- 
minosity functions of type 1 and type 2 quasars 
(Fig. [11]) and constrain the type 2 quasar fraction 
as a function of luminosity (Fig. [12|). We argue 
that the type 2 quasar luminosity function and the 
type 2 quasar fraction that we derive are robust 
lower limits. The main reasons are that (i) there 
may be a significant number of type 2 quasars that 
do not meet our spectroscopic selection criteria, 
and (ii) narrow lines are more extincted in type 2 
quasars than they are in type 1 quasars, biasing 
the type 2/type 1 ratio at a given luminosity to 
lower values. 

Objects with different [QUI] 5007 luminosities 
and different redshifts suffer from different se- 
lection biases. Our best data are at low red- 
shifts and relatively low luminosities (z < 0.3 and 
-^[Oiii] < 10^-^0) Eind at high redshifts and rel- 
atively high luminosities (0.50 < z < 0.83 and 
-^[oiii] > W^'^Lq). In these regimes we find that 
type 2 quasars are more abundant than type 1 
quasars, with the type 2/type 1 ratios of 1.5:1 and 
1.2:1, correspondingly. 
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Appendix 

In this appendix, we prove our assertion in 
H3.2.4I that the statistic u, defined by Eq. M fol- 
lows a distribution with one degree of free- 
dom. We have N quasars with the same redshift 
and [OIII]5007 luminosity, of which n are detected 
in the radio, so the measured detection rate is 
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/ = n/N . We would like to test the null hypothe- 
sis (NH) that the underlying probability of detec- 
tion (for this redshift and luminosity) is p (denoted 
as pkd{z, i[oiii]) in the text). If it were indeed p, 
then the probability of detecting n objects out of 
N would be binomial: 



i?(n) = C>"(l-p) 



N-n 



(16) 



For very large N {\/Np 3> 1), this distribution is 
close to a Gaussian 



r{n)dn 



1 



"dn, 



(17) 



where — Np(l — p). We now define a function 
of the observed values n and N and the underlying 
probability p: 



(n - pN)^ 
Np{l-p)' 



(18) 



The probability density of u, i.e., the probability 
to find this variable in the range between u and 
u + du under the assumption of the NH is 



ru{u)du 



2r{n) 
|du/dn 



■du 



1 



'2'KU 



-"/'dw (19) 



for w > 0, and r^iu) = otherwise. In other 
words, the distribution of u does not depend on 
N 01 p for large values of N. The distribution 
function ([T9l) is that of with one degree of free- 
dom. Given a measurement n and N , if we want 
to test whether it is consistent with the underlying 
probability p, we calculate u. If u > 6.5, then the 
NH is ruled out with a 99% probability. 

REFERENCES 

Abazajian, K., et al. 2003, AJ, 126,2081 (DRl pa- 
per) 

Abazajian, K., et al. 2004, AJ, 128, 502 (DR2 pa- 
per) 

Adams, T. F., & Weedman, D. W. 1975, ApJ, 199, 
19 

Adelman-McCarthy, J., et al. 2006, ApJS, 162, 1 
(DR4 paper) 

Adelman-McCarthy, J., et al. 2007, ApJS, 172, 634 
(DR5 paper) 



Adelman-McCarthy, J., et al. 2008, ApJS, 175, 297 
(DR6 paper) 

Anderson, K. 1970, ApJ, 162, 743 

Anderson, A., et al. 2007, AJ, 133, 313 

Antonucci, R. 1993, ARA&A, 31, 473 

Baldwin, J. 1975, ApJ, 201, 26 

Baldwin, J. A., PhiUips, M. M., & Terlevich, R. 
1981, PASP, 93, 5 

Barger, A. J., Cowie, L. L., Mushotzky, R. F., 
Yang, Y., Wang, W.-H., Steffen, A. T., Capak, 
P. 2005, AJ, 129, 578 

Becker, R. H., White, R. L., & Helfand, D. J. 1995, 
ApJ, 450, 559 

Beckmann, V., Soldi, S., Schrader, C. R., Gehrels, 
N., & Produit, N. 2006, ApJ, 652, 126 

Blanton, M. R., Lupton, R. H., Maley, F. M., 
Young, N., Zehavi, I., & Loveday, J. 2003, AJ, 
125, 2276 

Boroson, T. A. & Green, R. F. 1992, ApJS, 80, 
109 

Capriotti, E.R. 1964a, ApJ, 140, 632 

Capriotti, E.R. 1964b, ApJ, 139, 225 

Comastri, A., Setti, G., Zamorani, G. & Hasinger, 
G. 1995, A&A, 296, 1 

Dahari, O., & De Robertis, M. M. 1988, ApJS, 67, 
249 

De Zotti, G. & Gaskell, C. M. 1985, Astr. Ap., 
147, 1 

Eisenstein, D.J., et al. 2001, AJ, 122, 2267 

Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., 
Shimasaku, K., & Schneider, D. P. 1996, AJ, 
111, 1748 

GiUi, R., Comastri, A., & Hasinger, R. 2007, 
A&A, 463, 79 

Greene, J. & Ho, L. 2007, ApJ, 667, 131 

Grimes, J. A., Rawlings, S. & Willott, C. 2004, 
MNRAS, 349, 503 



19 



Groves, B., Dopita, M. & Sutherland, R. 2004, 

ApJS, 153, 75 

Gunn, J. E., et al. 1998, AJ, 116, 3040 

Gunn, J. E., et al. 2006, AJ, 131, 2332 

Haas, M., Siebenmorgen, R., Schulz, B., Krugel, 
E., & Chini, R. 2005, A&A, 442, 39 

Hao, L., et al. 2005a, AJ, 129, 1783 

Hao, L., et al. 2005b, AJ, 129, 1795 

Hasinger, G. 2008, A&A, 490, 605 

Heckman, T. M., van Brcugcl, W. J. M., & 
Butcher, H. R. 1981, ApJ, 247, 403 

Heckman, T. M., et al. 2005, ApJ, 634, 161 

Hopkins, P. F., Hcrnquist, L., Cox, T. J., Di Mat- 
teo, T., Robertson, B., & Springel, V. 2006, 
ApJS, 163, 1 

Hopkins, P. F., Richards, G. T., & Hernquist, L. 
2007, ApJ, 654, 731 

Huchra, J. & Burg, R. 1992, ApJ, 393, 90 

Ivezic, Z., et al. 2004, Astron. Nachr., 325,583 

Kauffmann, G., et al. 2003, MNRAS, 346, 1055 

Kewley, L. J., Dopita, M. A., Sutherland, R. S., 
Heisler, C. A., & Trevena, J. 2001, ApJ, 556, 
121 

Krolik, J.H. 1999, Active Galactic Nuclei: From 
the Central Black Hole to the Galactic Environ- 
ment (Princeton: Princeton University Press) 
[pp. 456-459 for biases in the line-selected sam- 
ples] 

Krolik, J.H. & McKee, C.F. 1978, ApJS, 37, 459 

Lacy, M., Canalizo, G., Rawlings, S., Sajina, A., 
Storrie-Lombardi, L., Armus, L., Marleau, F. 
R., & Muzzin, A. 2005, Memorie della Societa 
Astronomica Italiana, 76, 154 

Lawrence, A. 1991, MNRAS, 252, 586 

Lupton, R.H., Gunn, J.E., & Szalay, A.S. 1999, 
AJ, 118, 1406 



Marconi, A., Risaliti, G., Gilli, R., Hunt, L.K., 
Maiolino, R. & Salvati, M. 2004, MNRAS, 351, 
169 

Markwardt, C. B., Tueller, J., Skinner, G. K., 
Gehrels, N., Barthelmy, S. D., & Mushotzky, 
R. F. 2005, ApJ, 633, 77 

Martmez-Sansigre, A., et al. 2006, MNRAS, 370, 

1479 

Martmez-Sansigre, A., et al. 2007, MNRAS, 379, 
L6 

Mulchaey, J.S., Koratkar, A., Ward, M.J., Wilson, 
A.S., Whittle, M., Antonucci, R. R. J., Kinney, 
A. L., & Hurt, T. 1994, ApJ, 436, 586 

Nagao, T., Maolino, R. & Marconi, A. 2006, A&A, 

447, 863 

Netzer, H. 1975, MNRAS, 171, 395 

Netzer, H., Mainieri, V., Rosati, P., & Trakhten- 
brot, B. 2006, A&A, 453, 525 

Osterbrock, D. E. 1977, ApJ, 215, 733 

Osterbrock, D. E. 1989, Astrophysics of Gaseous 
Nebulae and Active Galactic Nuclei (Mill Valley 
CA: University Science Books) 

Osterbrock, D. E., Koski, A.T., & Phillips, M.M. 

1976, ApJ, 206, 898 

Osterbrock, D. E. & Shaw, R. A. 1988 ApJ, 327, 
89 

Padmanabhan, N., et al. 2008, ApJ, 674, 1217 

Peacock, J.A. 1983, MNRAS, 202, 615 

Petrosian, V. 1976. ApJ, 210, L53 

Pier, J. R., Munn, J. A., Hindsley, R. B., Hen- 
nessy, G. S., Kent, S. M., Lupton, R. H., & 
Ivezic, Z. 2003, AJ, 125, 1559 

PoUetta, M. 2006, ApJ, 642, 673 

Polletta, M., Weedman, D., Hoenig, S., Lonsdale, 
C. J., Smith, H. E., & Houck, J. 2008, ApJ, 
675, 960 

Ptak, A., Zakamska, N. L., Strauss, M. A., Krolik, 
J. H., Heckman, T. M., Schneider, D. P., & 
Brinkmann, J. 2006, ApJ, 637, 147 



20 



Richards, G. T., et al. 2001, AJ, 21, 2308 

Richards, G. T., et al. 2002, AJ, 123, 2945 

Richards, G.T., et al. 2006, AJ, 131, 2766 

Rigby, J. R., Rieke, G. H., Donley, J. L., Alonso- 
Herrero, A., & Perez-Gonzalez, P. G. 2006, 
ApJ, 645, 115 

Risaliti, G., Maiohno, R., & Salvati, M. 1999, ApJ, 
522, 157 

Salzer, J. 1989, ApJ, 347, 152 

Sandage, A. 1961, ApJ, 133, 335 

Sazonov, S., Revnivtsev, M., Krivonos, R., Chu- 
razov, E., & Sunyaev, R. 2007, A&A, 462, 57 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 
1998, ApJ, 500, 525 

Schmidt, M. 1968, ApJ, 151, 393 

Schneider, D. P., et. al. 2007, AJ, 134, 102 

Shankar, F., Weinberg, D. H., & Miralda-Escude, 
J. 2007, ApJ, submitted (|arXiv:0710.4488p 

Shen, Y., et al. 2007, AJ, 133, 2222 

Simpson, C. 1998, MNRAS, 297, L39 

Simpson, C. 2005, MNRAS, 360, 565 

Smith, J. A., et al. 2002, AJ, 123, 2121 

Soltan, A. 1982, MNRAS, 200, 115 

Stern, D., et al. 2005, ApJ, 631, 163 

Stoughton, C., et al. 2002, AJ, 123, 485 

Strauss, M. A., et al. 2002, AJ, 124, 1810 

Sturm, E., et al. 2006, ApJ, 642, 81 

Subbarao, M., Frieman, J., Bernardi, M., Loveday, 
J., Nichol, B., Castander, F., & Meiksin, A. 
2002, SPIE, 4847, 452 

Szokoly, G. P., et al. 2004, ApJS, 155, 271 

Tucker, D. L., et al. 2006, Astron. Nachr. 327, 821 

Treister, E., Urry, C. M., & Lira, P. 2006, Revista 
Mexicana de Astronomia y Astrofisica, 26, 147 



Treister, E., Krolik, J., & Dullemond, C. 2008, 
ApJ, 679, 140 

Ueda, Y., Akiyama, M., Ohta, K., & Miyaji, T. 
2003, ApJ, 598, 886 

Vanden Berk, D.E., et al. 2001, AJ, 122, 549 

Veilleux, S., & Osterbrock, D. E. 1987, ApJS, 63, 
295 

Vignah, C., Alexander, D., & Comastri, A. 2004, 
MNRAS, 354, 720 

Vignali, C., Alexander, D., & Comastri, A. 2006, 
MNRAS, 373, 321 

Voges, W., et al. 1999, A&A, 349, 389 

Urry, C. M. & Padovani, P. 1995, PASP, 107, 803 

Whittle, M. 1985, MNRAS, 216, 817 

York, D. G., et al. 2000, AJ, 120, 1579 

Yu, Q. & Tremaine, S. 2002, MNRAS, 335, 965 

Zakamska, N. L., et al. 2003, AJ, 126, 2125 (Paper 
I) 

Zakamska, N. L., Strauss, M.A., Heckman, T.M., 
Ivezic, Z., & Krolik, J.H. 2004, AJ, 128, 1002 

Zakamska, N. L., et al. 2005, AJ, 129, 1212 

Zakamska, N. L., et al. 2006, AJ, 637, 147 

Zakamska, N. L., et al. 2008, submitted 



This 2-column preprint was prepared with the AAS lAT^rjX 
macros v5.2. 



21 



Table 1 

Catalog of 887 Optically-Selected Type 2 Quasars 



RA 




Dec 


redshift 


,(G) 
[OIIIl 


(NP) 
[OIIIl 


Plate 


Fiber 


MJD 


Tcode 


S20cm 


RMS 


"PSF 


3PSF 


'■PSF 


•PSF 


^PSF 




"a 








0.746236 





671701 


0.6007 


9.04 


9.12 


686 


350 


52519 


00001 


0.00 


0.149 


21.67 


21.35 


20.54 


19.72 


19.83 


0.36 


0.06 


0.04 


0.03 


0.11 


0.965958 


-1 


028334 


0.2689 


8.16 


8.15 


669 


289 


52559 


00001 


8.45 


0.148 


21.93 


20.68 


19.26 


18.84 


18.31 


0.20 


0.03 


0.02 


0.02 


0.03 


1.222618 


-0. 


843989 


0.6430 


9.18 


9.20 


669 


209 


52559 


00001 


0.00 


0.153 


22.37 


21.90 


21.03 


20.19 


19.97 


0.32 


0.11 


0.06 


0.04 


0.13 


1.870121 


1 


,101123 


0.4663 


8.33 


8.33 


669 


457 


52559 


00001 


0.00 


0.084 


21.80 


21.80 


20.29 


19.66 


19.19 


0.21 


0.07 


0.03 


0.03 


0.07 


2.799806 


0.940632 


0.4094 


8.67 


8.66 


669 


602 


52559 


00001 


0.00 


0.117 


21.89 


21.10 


19.84 


19.19 


18.99 


0.30 


0.06 


0.03 


0.02 


0.06 


2.862281 


15 


891562 


0.0999 


8.16 


8.17 


752 


380 


52251 


10000 


0.00 


-1.000 


20.72 


19.28 


18.50 


17.88 


17.61 


0.09 


0.02 


0.02 


0.02 


0.03 


3.026290 


-9. 


790457 


0.1668 


8.54 


8.55 


652 


399 


52138 


10000 


0.00 


0.153 


20.42 


18.95 


18.16 


17.81 


17.69 


0.07 


0.02 


0.02 


0.02 


0.03 


6.070319 


-9. 


545760 


0.3600 


8.60 


8.62 


1913 


381 


53321 


00001 


0.00 


0.149 


20.05 


19.46 


18.74 


18.54 


18.27 


0.05 


0.02 


0.01 


0.02 


0.03 


5.344129 


-0. 


254834 


0.5493 


8.35 


8.36 


687 


22 


52518 


00001 


1.48 


0.175 


21.35 


21.02 


20.24 


19.82 


19.61 


0.15 


0.04 


0.03 


0.03 


0.08 


6.381066 


-10. 


672835 


0.3035 


8.64 


8.65 


653 


149 


52145 


01000 


1.16 


0.135 


19.78 


19.59 


19.06 


18.94 


18.50 


0.05 


0.02 


0.02 


0.02 


0.04 



Note.— (1) RA and Dec arc in .72000.0 coordinates. 

(2) ^Iqjjj] and are listed in units of log(i/LQ) and refer to Gaussian and non-parametric measures of the [OIII]5007 line luminosity, respectively 

(JHH. 

(3) The first four digits of the target code (col. 9) show whether the object was targeted with the Galaxy, Low-z QSO, High-z QSO, or Serendipity FIRST 
algorithm, in that order. The last digit indicates whether the object belongs to the Special Southern survey. 

(4) '520cm the peak flux (in ni.Jy/beam) of the nearest match within '2" from the FIRST survey jBecker et a1.|P^^ . and is listed as zero, if there is no 
match. If the object is not in the FIRST survey area, the RMS flux value is listed as —1. 

(5) Columns 12-16 and 17-21 list ugriz PSF asinh magnitudes from the SDSS TARGET database, corrected for Galactic extinction jgchlegel et al.|P^^ . 
and their 1-sigma errors. 

Table 1 is presented in its entirety in the electronic edition of the Astronomical Journal. A portion is presented here for guidance regarding its form and 
content. 



Table 2 

SDSS Spectroscopic Target Algorithms 



Target Algorithm 


No. 


of Obj. 


t^aiti Survey 


Main Galaxy 




366 


Main Low-z QSO 




76 


Main High-z QSO 




89 


Main Serendipity FIRST 




276 


Main LRG 




64 


Main QSO FIRST 




14 


Khun li( )SAT 










771 


Special Southern Survey 


Southern Galaxy 




4 


Southern LRG 




5 


Southern Low-z QSO 




20 


Southern High-z QSO 




6 


Southern Serendipity FIRST 




9 


Southern ROSAT 




1 


Faint quasars 




31 


Photo-z 




29 


Faint LRG 




5 


w-band galaxy 




2 


Oombined 




116 


Note.— SDSS spectroscopic 


target algo- 


rithms important for select 


,ing 


type 2 


quasars, and the number of obj 


ects 


targeted 


by each. Objects can be targi 




by mul- 


tiple algorithms, so these nur 


Qbcr 




add up to the total number of 


objc 


lets. The 


type 2 quasar luminosity functi 


on p 


ircscntcd 


in this paper is derived from th 


c 740 objects 


targeted by the top four Main 


target algo- 


rithms: Galaxy, Low-z QSO, 


High-z QSO 


and Serendipity FIRST. 







Table 3 

Mean V/Vmax values calculated from the type 1 and type 2 quasar samples. 





No. 


(V/Knax> 


Type 


1 Quasars 




0.00 < z< 0.30 


1020 


0.55 


0.30 < z< 0.50 


2802 


0.56 


0.50 <z< 0.83 


4181 


0.57 


Combined 


8003 


0.61 


Type 2 Quasars 


0.00 <z< 0.30 


420 


0.46 


0.30 <z< 0.50 


175 


0.61 


0.50 < z < 0.83 


145 


0.54 



Note. — Our type 1 quasar sample 
is complete and values greater than 0.5 
indicate that there is positive redshift 
evolution within each redshift range. 
Our type 2 quasar sample is incom- 
plete and the deviation of (V/Vma-x) 
from 0.5 is due to a combination of 
redshift evolution and selection effects, 
which we do not attempt to disentan- 
gle. 
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4900 5000 6500 6600 6500 6600 
rest-frame wavelength, A 

Fig. 1. — Examples of our non-parametric fitting 
procedure to determine whether there is a broad 
component in Ha hue. From left to right, the three 
panels for each object are: (i) the H/3-|-[0III] com- 
plex and model fits to the blue wing of [OIII]4959 
and the red wing of [OIII]5007; (ii) fit of the 
Ha-h[NII] complex, using the [OIII]4959, 5007 hne 
profile for all three lines; (iii) the best 4-Gaussian 
fit to the Ha-|-[NII] complex. In each panel, the 
original spectrum is in black, the model is in red, 
residuals are in blue and the four Gaussian com- 
ponents are in green. Although for both objects 
the 4-Gaussian fit to the lla+ [Nil] complex is sta- 
tistically preferred to the 3-Gaussian fit, in SDSS 
JllOO+0846 (upper panel) this complex is well-fit 
with a blend of lines shaped like [OIII] , so we keep 
this object in our sample. In SDSS J1648-h3022 
(lower panel), the complex is significantly broader 
than the blend of [OIII] lines, as indicated by the 
fit residuals, so we classify this object as a broad- 
line AGN and exclude it from our sample. 



Fig. 2. — Example spectra of the highest luminos- 
ity type 2 quasars in the redshift range z < 0.30, 
0.30 < z < 0.50, and 0.50 < z < 0.83 smoothed 
by 5 pixels for display purposes. Objects in the 
three panels have i[oiii]> 1-8 x 10^, 5.0 x 10^, and 
9.6 X 10^ Lq, respectively. The strongest emission 
line is [OIII]5007. 
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0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 
redshift 

Fig. 3. — Distribution of [OIII]5007 luminosities 
and redshifts of type 2 quasars targeted by the 
Galaxy (366 objects), Low-z QSO (76), High-z 
QSO (89), Serendipity FIRST (276), and other 
(147) target algorithms. Some objects are tar- 
geted by multiple algorithms and thus appear in 
multiple panels in this Figure. Also shown is the 
distribution for type 1 quasars (8003 objects; see 
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Fig. 4. — Comparison of the observed probabil- 
ity of FIRST detections with the best-fitting func- 
tion for Prd(^, i[oiii]) (Eq- rroi) for non-radio- 
selected type 2 quasars. Left: in each redshift 
range zi < z < Z2, the circle shows the value 
n/N, where n is the number of radio detected ob- 
jects in this bin and N is the total number of ob- 
jects. The Poisson error bars reflect the total num- 
ber of objects contributing to each bin. The solid 
line connects values Pave{z) = '}2iPm{zi,Li)/N , 
where the summation is over all objects in this 
bin. Right: same, but for luminosity bins instead 
of redshift bins. Such representation allows us to 
properly take into account the distribution of ob- 
jects in the z-L^qyo] plane. 




log(L[oiii]/Lo) 

Fig. 5. — [OIII]5007 luminosity function of type 
2 quasars (lower limits) for three redshift ranges, 
with and without correction for the probability 
of radio detection (see H3.2.41 solid and dashed 
curves, respectively). The maximum volume is 
calculated using the selection criteria of the var- 
ious SDSS spectroscopic target selection algo- 
rithms (see t j3.2p . Included in the calculation are 
871 objects targeted by the top four Main survey 
algorithms (Galaxy, Low-z QSO, High-z QSO, and 
Serendipity FIRST). 
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O Hao et al. (2005b) 
• This work 



10= 



Fig. 6. — [OIIIJ5007 luminosity function of type 
2 quasars (same as blue solid curve in Fig. [5]) 
for z < .30, co mpared with that derived by 
(l2005bl) . 



Hao et al 



We find good agreement be- 
tween the two functions in the luminosity range in 
which they overlap. Here, the luminosities from 
Hao et al.l (|2005bl had been shifted up by 0.14 dex 
to account for the difference in their spectropho- 
tometric flux calibration scale (see 
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Fig. 7. — Contributory power of the various SDSS 
target algorithms (as labeled)- the fraction of ob- 
jects in the type 2 quasar sample that would be se- 
lected by the algorithm taken alone. We place ob- 
jects at different redshifts (in steps of Az — 0.01) 
to determine this quantity as a function of as- 
sumed redshift. The three panels show the results 
for different ranges in [OHIjSOOT luminosity. All 
four algorithms are poor at selecting low i[oiii] 
objects at high redshifts, which is why the derived 
type 2 quasar LF shows evidence for incomplete- 
ness there. 



26 




l°8(L[oiii]/Lo) 

Fig. 8. — [OIII]5007 luminosity function of type 
1 quasars for three ranges in redshift. The max- 
imum volume is calculated using the magnitude 
limit of the Low-z QSO target algorithm. The 
turn-over of the 0.50 < z < 0.83 LF at low lu- 
minosities is an artifact due to the difficulty of 
measuring weak [OIII]5007 lines at these high red- 
shifts. 




Fig. 9. — Correlation between broadband abso- 
lute magnitude M2500 and [OIII]5007 luminosity 
for type 1 quasars in our sample. M2500 measures 
the continuum luminosity at around 2500 A, cor- 
responding to the SDSS i band filter at z — 2. We 
find the mean relation to be log(L[OIII]/Lo) = 
— 0.38Af25oo— 0.62 (dashed line), i.e., the relation is 
close to linear. The scatter in i[oiii] at fixed con- 
tinuum luminosity is consistent with a log-normal 
scatter with a width of 0.36 dex; 1-sigma and 2- 
sigma contours are shown here (solid curves). 
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Fig. 10.^ Comparison of the derived [OIII]5007 
luminosity function of type 1 quasars with pre- 
vious wo rk. We convolv e the broadband LF de- 
rived by [Richards et al. 



(120061 ) for 0.30 < 



z < 



0.68 with the mean i[oiii]--^2500 relation Eq. [12] 
with a log-normal scatter of 0.36 dex. The 
converted [OIII]5007 luminosity function (short- 
dashed curve) is in good agreement with the de- 
rived LF in the redshift range 0.30 < z < 0.50. 
The shape of the converted LF is sensitive to the 
assumed scatter in the i[oiii]--^25oo relation. Also 
shown is the converted [OIII]5007 LF derived un- 
der the assumption of zero scatter (long-dashed 
curve), which is inconsistent with the derived LF. 
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Fig. 11. — Comparison of the [OIIIJ5007 luminos- 
ity functions of type 1 and type 2 quasars, for 
three ranges in redshift. The type 2 quasar lumi- 
nosity function (lovi^er limits) with and without 
the correction for the probability of radio detec- 
tion (see §3.2.4|) are shown (black solid and blue 
open circles, respectively). The type 2 quasar sam- 
ple suffers from incompleteness, especially at low 
[OIII]5007 luminosities and high redshifts. Nev- 
ertheless, we find that the derived space densities 
of type 1 and type 2 quasars are comparable for 
the redshift range z < 0.30 and for the highest lu- 
minosities for the higher redshifts 0.3 < z < 0.83. 
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Fig. 12. — Lower limits to the type 2 quasar frac- 
tion: the ratio of type 2-to total (type 1 + type 2) 
quasar number densities for three ranges in red- 
shift: z < 0.30 (blue), 0.30 < z < 0.50 (green), 
and 0.50 < z < 0.83 (red). Number densities 
are estimated by integrating the LF over bins in 
[OIII]5007 luminosity. We have used our best 
lower bound to the type 2 quasar LF (which in- 
cludes a correction for the probability of radio de- 
tection; see Section 3.2.4 in the main paper). Ob- 
scured quasar fractions derived from X-ray surveys 
are from Ueda et al. (2003; grey open squares) 
and Hasinger (2008; black open circles and dashed 
line); those derived from IR data are Treister et al. 
(2008; black filled triangles). Also shown are re- 
sults from radio data (Grimes et al. 2004; grey 
open triangles) and optical data (Simpson et al. 
2005; grey open diamonds). 
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Fig. 13. — Balmer line ratios of luminous type 2 
quasars, for objects in which the relevant emission 
lines can be measured. Left panel: H/3/H7 vs. 
[OIII]5007 luminosity (203 objects). The median 
H/?/H7 is 3.06 (dotted line), which corresponds to 
an [OIII]5007 extinction of 2.6 mag (for a Milky 
Way extinction curve and a foreground obscuring 
screen). Right panel: Ha/H/3 vs. H/3/H7 (44 ob- 
jects). The arrow shows the extinction vector for 
the Milky way extinction curve, whose length cor- 
responds to an [OIII]5007 extinction of 2.2 mag. 
In both panels, the hatched area shows case A-B 
ratios at temperatures 5,000-20,000 K, and errors 
are estimated by comparing non-parametric and 
Gaussian flux measurements. 
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